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ABSTRACT 


The conventional linear distribution of stress in a section 
subjected to bending in the plastic range, as usually represented 
by the classical my/I formula, has been replaced by a simple 
trapezoidal distribution that varies with the applied bending 
moment. This latter distribution approximates the distribution 
based on the actual stress-strain characteristics of the material 
and the shape of the cross section. 

A method has been developed for use throughout the elastic 
and plastic range of bending which covers the following points: 
(1) the resisting moment of any cross section corresponding to 
any extreme fiber stress up to the ultimate stress of the material; 
(2) the corresponding stress distribution over the entire cross 
section; (3) the corresponding shear distribution; (4) the effects 
of local crippling; (5) complex bending; and (6) the combined 
effects of axial load and bending. 


SyMBOLsS USED 


Main (as distinguished from subscript) 


A = area between fiber in question and the extreme fiber 
a = angle the neutral axis makes with respect to the 
X axis 


= maximum width of the section 


B = angle the moment vector makes with respect to 
the X axis 
c = perpendicular distance from a principal axis to the 


extreme fiber 
e = unit strain 
f = applied stress 
F = allowable stress 
¢ = angle the principal X’ axis makes with the X axis 
7 = an additional (multiplying) factor of safety 
k = section factor 
m = applied moment 


M = allowable moment 
r = normal load on portion of cross section 
Q = static (first) moment, about the neutral axis, of 


the area between the fiber in question and the 
extreme fiber 

R = rate of change of trapezoidal stress intercept with 
respect to the stress of the extreme fiber; or 
stress ratio 

S$ = shear load 
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= distances measured parallel to the Y and Y axes 


x,y 

X, Y = mutually perpendicular axes 

X’, Y’ = principal axes 

Z = width of section at distance y from the neutral axis 

Subscript 

a = due to axial or direct load 

b = due to conventional bending 

d = for a diamond type of cross section 

h = for an hour-glass type of cross section 

L = due to limit or applied load 

m = maximum 

n = due to normal or conventional practice 

0 = minimum 

r = for a rectangular type of cross section; or resultant 

t = tensile or total 

u = ultimate 

x = about X axis 

y = about Y axis; or yield; or fiber at distance y from the 
neutral axis 

x’ = about X’ axis 

y’ = about Y’ axis 


INTRODUCTION 


hen CLASSIC BENDING THEORY is based on the 
assumption that plane sections before bending 
remain plane after bending. Since, within the elastic 
range, stress is proportional to strain, it follows that the 
stress distribution over the cross section will be linear, 
varying from zero at the neutral axis to a maximum 
at the extreme fiber. This distribution is expressed 
mathematically for a symmetric cross section subjected 
to simple bending by the conventional formula 


fo = my/I (1) 
where 


fy = bending stress at a fiber that is located at a 


distance y from the neutral axis 
the applied moment (lower case used to differ- 


m = 
entiate applied moment m from allowable 
moment M) 

I = the moment of inertia about the axis of bend- 
ing 
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When, on the other hand, the extreme fiber stress 
exceeds the proportional limit stress of the material, 
the stress distribution over the cross section is no longer 
linear but tends to follow the stress-strain characteris- 
tics of the material.! Consequently, the conventional 
formula for linear bending stress distribution no longer 
applies. It has been standard aircraft practice, in 
such cases, to use a “form factor’ to indicate the 
amount by which the allowable ultimate stress of the 
material could be considered to be increased in 
order to apply the conventional f, = my/J form- 
ula. 

The use of this simple expedient of a form factor 
necessitates having available the form factors for all 
the various shapes of cross sections commonly en- 
countered in practice. Such form factors for round 
tubes have been incorporated in the “Bending Modulus 
of Rupture Curves for Round Chrome Molybdenum 
Steel Tubing.’ Form factors for sections other than 
round are not available to the average structural 
engineer. In view of this lack of material, it was 
proposed to make a further study of form factors for 
various I-type, channel-type, and round tubular types 
of cross sections. Simple bending tests were made 
on such sections, and the data’ pertinent to this article 
are included in the Appendix. 

The current ‘‘form factor’’ method of accounting for 
the effects of bending in the plastic range is subject 
to the following limitations: 

(1) The physical significance of the plastic bending 
phenomenon is lost by the use of a fictitious stress 
greater than the allowable ultimate stress F,, of the 
material. 

(2) In many cases the full advantage of the form 
factor effect at ultimate loads cannot be realized since 
the yield strength criteria must also be met, particularly 
in aircraft work. Therefore, yield form factors as 
well as ultimate form factors should be made avail- 
able. 

(3) The stress distribution across the entire cross 
section is fictitious and, therefore, local crippling 
effects cannot be taken care of unless a sufficient number 
of tests are run on various sections of varying wall 
thickness. (The bending modulus of rupture curves 
for round tubes of various D/t values include this 
“local crippling” effect.*) 

(4) In order to include the effect of material proper- 
ties, a complete set of form factors would have to be 
made available for each material. 

To remove the above limitations it was decided to 
abandon the conventional form factor theory and 
evolve another method for handling the plastic bending 
phenomenon. The use of a bending stress distribution 
based on the exact stress-strain characteristics of the 
material was considered but was found to be imprac- 
tical because of the difficulty of expressing the stress- 
strain relationship mathematically. Beilschmidt,* how- 
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ever, has approximated this relationship by expressing 
stress f as a function of some power of strain e, (f = 
ae"). This method, however, was also limited by the 
mathematical difficulties encountered in applying it 
to typical problems. 

With the idea of utilizing, whenever possible, the 
common physical properties of a cross section, a 
trapezoidal distribution approximating the actual 
stress-strain distribution was found to give a practical 
solution. The trapezoidal distribution was varied, of 
course, depending on the extreme fiber stress in- 
volved. 

A formula was developed to determine the moment 
of resistance offered by a cross section subjected to 
any bending moment m. The corresponding formula 
for shear distribution was also developed. 

In order to account for the effect of material proper- 
ties, curves enabling the designer to account for plastic 
bending phenomena were developed for each of various 
materials commonly used in aircraft practice. 

A discussion of the proposed method of handling 
local crippling, complex bending, and combined axial 
load and bending is included, together with the re- 
spective limitations of the proposed method. 


MOMENT OF RESISTANCE OF A SECTION SUBJECTED 
TO BENDING 


Linear or Triangular Stress Distribution 


The stress-strain curve may be thought of as a stress 
distribution curve if the vertical axis is assumed to be 
the neutral axis. 

Fig. 1(a) shows a typical stress-strain curve in which 
the strain of the extreme fiber e, corresponds to a stress 
fy which is below the proportional limit of the material. 
The equivalent stress distribution is shown on section 
A-A of Fig. 1(b). 


WA 
Stress, 
* 








Fic. 1. Triangular stress distribution. 


The moment of resistance m, of the section subjected 
to bending may be expressed as 


m, = fr(I/c) (2) 


Subscript 6 designates a triangular distribution; 
fy is the corresponding extreme fiber stress; and I/c 
is the conventional section modulus. 
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Constant or Rectangular Stress Distribution 


If a stress-strain curve of a material were as shown in 
Fig. 2(a), the shaded area could be thought of as the 
stress distribution. The equivalent stress distribution 
is shown on section A-A of Fig. 2(b). 

The moment of resistance m, of the section subjected 
to bending may be expressed approximately as 


My = 2f:Om (3) 


Subscript r designates a rectangular distribution; 
f, is the corresponding extreme fiber stress; and Q,, is 
the maximum static moment of area (i.e., the first 
moment about the neutral axis, of the entire area 
between the neutral axis and the extreme fiber). 





Fic. 2. Rectangular stress distribution. 


Trapezoidal Stress Distribution 


Consider the condition where the extreme fiber stress 
fm of a section subjected to bending lies between the 
allowable yield stress F,, and the allowable ultimate 
stress F,,, as shown in Fig. 3(a). The stress-strain 
distribution has been approximated by a trapezoidal 
distribution as indicated by the shaded area. The 
stress indicated by the intercept of the trapezoid with 
the vertical axis through the origin is designated as the 
trapezoidal intercept stress fy. The equivalent trape- 
zoidal distribution is shown on section A-A of Fig. 
3(b). 








Fic. 3. Trapezoidal stress distribution. 
The moment of resistance m of the section subjected 
to bending may be expressed as 
m= M + m, (4) 


Since the trapezoidal distribution may be considered 
as made up of a rectangular distribution with a maxi- 


mum fiber stress of f, and a triangular distribution with 
a maximum fiber stress of f). 

The values of m, and m, when substituted in Eq. 
(4) will give 


m= fi(T/c) + 2f,0m (5) 


However, it may be seen from Fig. 3(a) that f, = fo 
the trapezoidal intercept stress and that f, = fm ~ fo. 

Therefore, these values when substituted in Eq. (5) 
will give 


m= (Fo —_ fo Cc) + 2foO m (6) 


By regrouping the terms, this expression may be 
further reduced to 


m = fn(I/c) + fol2Qn — (I/c)] (7) 


In order to utilize the conventional expression mc/I, 
divide through by //c to obtain 


me _ 20m _ 
Ba +1(7 ) = 

The term mc/I may be thought of as a measure of 
the moment of resistance that could be developed by 
the section, for when multiplied by the section modulus 
I/c it will give the moment of resistance. 

The term 2Q,,/(J/c) which depends only on the shape 
of the cross section is designated by k and will be re- 
ferred to as the “‘section factor.” 

In its final form the expression reduces to 


mc/I = fm + folk — 1) (9) 


Basic Assumptions 


In order to apply the expression for trapezoidal dis- 
tribution to the general problem of plastic bending, the 
following assumptions are made: 

(1) The tensile and compressive stress-strain proper- 
ties of the material are identical. This gives a simple 
distribution and avoids the difficulty of a complex 
distribution. These difficulties will be discussed in 
the latter part of this article. 

(2) The trapezoidal intercept stress fy is to be 
determined so that the first moments about the neutral 
axis of the area under te actual stress distribution 
curve are equal to the first moments about the neutral 
axis of the area under the equivalent trapezoidal stress 
distribution curve. It is true that such an approxi- 
mation is exact for a rectangular cross section only. 

The variety of cross sections encountered in practice 
precluded an investigation on all of them. Therefore, 
it was decided to consider that most aircraft sections 
could be broken down into three types. " 

Type (a) is the rectangular type, in which the 
material is distributed uniformly from the neutral axis 
to the extreme fiber. 

Type (b) is the diamond type in which the amount 
of material decreases from a maximum at the neutral 
axis to a minimum at the extreme fiber. Under this 
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classification would come the solid round and relatively 
thick-walled tubes. 

Type (c) is the more common hour-glass type, in 
which the amount of material increases from a mini- 
mum at the neutral axis to a maximum at the extreme 
fiber. Under this classification would come the J, 
channel, and most of the tubular sections. 

The relative values of trapezoidal stress intercept 
for these three types have been shown, in the Appendix, 
to be approximately equal. This implies that the 
resisting moments for all three types can be predicted 
by the use of the equivalent trapezoidal distribution 
based on a rectangular cross section. 














Fic. 4. Type (a). Rectangular type. 





Diamond type. 


Fic. 5. Type (b). 








Fic. 6. Type (c). Hour-glass type. 


Limits of Factors 


The range of values for the section factor k for various 
sections is shown in Table 1. 

It may be interesting to note that Table 1 will be 
found useful to check preliminary section properties. 
If the value for a section factor lies outside of the range 
shown in Table 1, it is fairly certain that an error has 
been made. This simple check has been used 
extensively in the Lockheed Engineering Depart- 
ment. 


SHEAR DISTRIBUTION 


The shear stress distribution corresponding to the 
trapezoidal bending stress distribution will show a 
marked change over that obtained from the con- 
ventional formula: , 


f, = SQ/dl (10) 
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TABLE 1 _—s 

Values of Section Factor, k = 2Q,,/(I/c) 
Section Value of k 
— t- 1.0 
— g- 1.5 
- Ee - 1.0 to 1.5 
- _— 1,0 to 1.5 





1.25 to 1.70 


_ @- 2.0 


1.5 to 2.0* 


— 4 = 4/3 


* Any section having a k value larger than 2.0 should be used 
with caution as it indicates a relatively thin outstanding leg. 























Since it is more convenient to work with shear 
flow, Eq. (10) will be replaced by the equation 


g = SQ/I 


where g equals the shear load per inch. 
The new expression for g is derived in the Appendix 
and may be written as 


_ SQfl + Ri(Ac/Q) — 1) | 
2 ‘ a 


where Q equals the first moment about the neutral 
axis, of the area A between the fiber in question and 
the extreme fiber and R equals df,/af, or the rate of 
change of the trapezoidal intercept stress with respect 
to the maximum stress. 

In the elastic range where the trapezoidal distribution 
becomes triangular (fp = 0), the value of R is zero and 
the expression for shear flow reduces to the conventional 


form. 


(10a) 


1+ R(k — 1) 


DEVELOPMENT OF CURVES OF PLASTIC BENDING AND 
SHEAR 


In order to make use of the relations discussed in the 
previous pages, certain factors and properties of 
individual materials will have to be made available. 

The basic need is for a complete stress-strain diagram 
for each material. By complete is meant a diagram 
that shows the stress-strain properties up to and 
including the ultimate stress. As the aeronautical 
structural engineer is faced with more and more re- 
finement of design, the need for complete stress-strain 
diagrams becomes more urgent. In fact, much of the 
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future development in structural analysis will be 
dependent on the plastic stress-strain characteristics 
of materials. 

Stress-strain curves have been assumed for some 
common aircraft materials so as to satisfy the minimum 
guaranteed values of: (1) modulus of elasticity, (2) 
proportional limit, (3) yield stress, (4) ultimate stress, 
and (5) ultimate elongation. 

By means of a tabular form, a copy of which is shown 
in the Appendix, the stress-strain curve for each of 
several materials was integrated numerically to give 
an equivalent trapezoidal distribution with a minimum 
stress at the neutral axis designated as fy and a maxi- 
mum stress at the extreme fiber designated as fy. 
Values of fm were assumed at various points on the 
stress-strain curve, and the corresponding fy was found 
so that the first moments about the vertical axis through 
the origin were equal for (1) the area under the trapezoid 
and (2) the area under the stress-strain curve. This 
value of f) was substituted into the expression: 


mc/I = fm + folk — 1) (12) 


for values of k = 1.0, 1.2, 1.4, 1.6, 1.8 and 2.0. 

The values of intercept stress fo were plotted versus 
fm and the slope of this curve (dfo/dfm) was in turn 
computed numerically to give values of R. 

Then for each material, there appears on one sheet: 
(1) a plot of stress f versus strain e, (2) a plot of stress 
f versus mc/I for various values of k, (3) a plot of 
stress f versus R, and (4) values of fp/F,, and 
fo/ Fu: 

All the available curves of this type have been 
included as Figs. 9 through 17 of this article in the form 
in which they were used in the Lockheed Engineering 
Department. 


GENERAL DISCUSSION OF CURVES 
Bending Stress Distributions 


If the complete stress-strain curve of a material 
were available, any value of extreme fiber stress (f,,) 
within the limits of the material could be assumed 
and a complete picture of the actual bending stress 
distribution over the remainder of the cross section 
could be had. On the assumption that plane cross 
sections remain plane after bending, the stress at any 
intermediate point (y) corresponding to the propor- 
tionate strain (e,) of that point relative to the maximum 
strain (e,) of the extreme fiber (i.e., e, = (y/c)ém) 
could be determined. 


Moments of Resistance Corresponding to an Extreme 
Fiber Stress (f,,,) 


Corresponding to any extreme fiber stress (f,) the 
moment of resistance m may be found by entering the 
stress-(mc/I) curve with f, and determining the 
mc/I value for the proper ‘section factor” k. 

Then m = (mc/I)(I/c). 












































TABLE 2° 
Comparison of Test and Predicted Moments 
2 4 16 17 
e | Ratio 
al a 2 | Actual |Average 
s - g | TO Pre-| X-Ray 
$| 8 E dicted Rating 
© > | Moment 
a| § 2 15 for 
a sie Casting 
31] 1.104 - 
80} 1.105 - 
33 | 1.079 - 
13 | 1.064 - 
# 14] 1.070 - 
bo | 35) 1.045] - 
£& | 37] 1.039 - 
es 1 1.054 - 
st 2] 1.028 - 
3 3] 1.020 - 
8 4/ 1.059] = 
b 4 5] 1.030 - 
a 6 | 1.040 - 
an 7) 1.013 = 
8] 8| 1.040] - 
8% | 9] 1.057] - 
£& | 10] 1.062 - 
° 11 | 1.050 - 
= 112] 1.013] - 
wo| *8 1.114 | Fair- 
© 
& re 49 | 1.135 | Poor 
4 t2 75 | 12225 - 
“+o; 78 | 1.113 - 
@ a0 - 
1 HE 
1 ed Bed ie 
ole o| 17] 1.125 ~- 
o os 
| 19 | 1.002 - 
3/& £| 39 | 1.336 | Fair 
aes 40 | 1.131 | Fair- 
m1 S| 41 | 1.440 | Fair 
29 1.078 - 
o| | 27] 1.087| - 
§|& %| 21] 0.965| - 
B15 5125] 1.095] - 
° fe . 
2 23 | 1.024] - 
“ 62 | 1.054 | Poor 
E|w $8] 60 | 1.080 | Poors 
F: ‘| 61 | 1.235 | Fair- 
= s 59 | 1,075] Poor+ 
63-A} 1.117 | Poor 
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When the extreme fiber stress is chosen as the 
allowable yield stress F,,, then the resisting moment 
will be the allowable yield moment, /,. 

When the extreme fiber stress is chosen as the 
allowable ultimate stress F,,, then the resisting mo- 
ment will be the allowable ultimate moment, /,,. 

For any applied moment m the corresponding 
extreme fiber stress and the stress distribution over 
the entire cross section may be found by entering the 
curve with mc/I and determining the extreme fiber 
stress f,, for the proper k value of the section. After 
the extreme fiber stress has been obtained, the stress 
distribution over the remainder of the section may be 
found for any intermediate point as described under 
the heading ‘‘Bending Stress Distributions.” 


Comparison of Test Results 


Simple bending tests* had been run at Lockheed on 
various /, channel, and round sections for 14ST forging, 
magnesium forging, and 195-T6 castings. Ratios of 
actual test values to values computed by the proposed 
method are given in Table 2. A more complete table, 
showing the computations, is included in the Appendix. 

To reduce test results from actual physical properties 
to minimum guaranteed properties, a correction factor 
based on the relative first moments of the respective 
trapezoidal distributions, about the neutral axis, had 
to be determined. The derivation of this correction 
factor is also included in the Appendix. 

Correction factors were applied to all but the casting 
tests, since, in the case of casting design, additional 
factors of safety are included to take care of possible 
defects in castings. Physical properties of castings 
determined from a test billet do not necessarily predict 
the physical properties of the casting proper; therefore, 
average X-ray ratings for the castings are given in 
Column 17, Table 2. 

Column 16, Table 2, shows the ratio of actual 
allowable test moments (reduced to minimum guaran- 
teed values except for the castings) to allowable com- 
puted moments based on minimum guaranteed proper- 
ties. There is only one case in which this ratio is less 
than unity. Specimen 21, which is a channel section 
of 14ST forged material, shows a ratio of 0.965, which 
is equivalent to a negative margin of 3.5 per cent. 
Other specimens from the same billet show positive 
margins. 

The average ratio for the various type sections and 
materials is given in Table 3. g 


Margins of Safety 


In the plastic range, the margins of safety may no 
longer be found from the ratio of the allowable stress 
to the actual stress, since actual stress is no longer 
directly proportional to loading. The margin of safety 
must of necessity be found from the ratio of allowable 
moment to actual moment. 
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TABLE 3 
Average 
Section Material Ratio 
I 14ST forging 1.057 
Mg. forging 1.039 
195-T6 casting 1.147 
Round 14ST forging 1.078 
195-T6 casting 1.302 
Channel 14ST forging 1.040 
195-T6 casting 1.112 
Shear Flow 


The shear flow (g) for any section may be found by 
entering the curve with the mc/I value and determining 
the corresponding value of R for the proper k value of 
the section. This value of R when substituted in the 
expression for shear flow will give the desired result. 


Shift in Neutral Axis 


The actual neutral axis of the section will not always 
pass through the geometric centroid of the section. For 
instance, consider the case in which, although the 
tensile and compressive properties of the material are 
identical, the section is not symmetric about its cen- 
troidal axis. For such a case the neutral axis will no 
longer coincide with the centroidal axis in the plastic 
bending range. Although this effect is usually neg- 
ligible, the designer should be aware of its presence. 


Local Crippling 


A discussion of modulus of rupture effects would not 
be complete without mention of local crippling phenom- 
ena. This is especially true in connection with 
structural aircraft practice, since the “Bending Modulus 
of Rupture Curves for Round Chrome Molybdenum 
Steel Tubing’? has been in use for many years. 

An application of the proposed method to ‘“‘chrome- 
moly”’ steel tubes heat-treated to 150,000 Ibs. per sq.in., 
for various ratios of D/T gave results shown by curve 
AFBC in Fig. 7. Curve DBGE has been reproduced 
from Figs. 4-20 of ANC-5.’ 

On the surface it would seem that the two methods 
are inconsistent. The differences, however, in the 
ranges of D/t larger than 12 may be attributed to the 
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Fic. 7. Bending modulus of rupture ‘‘chrome-moly”’ steel heat- 


treated to 150,000 Ibs. per sq.in. 
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effects of local crippling. There is a definite need, 
therefore, for actual crippling data which, if available, 
could be used as allowable fiber stresses in the proposed 
method to predict allowable bending moments. For 
the present it will be conservative to use the proposed 
method and ensure that no fiber will be stressed beyond 
its ‘‘crippling stress in pure compression,’’ whenever 
such data are available. For round tubes, however, 
the ANC-5 bending modulus of rupture curves should 
be used, since they already include the crippling 
effects. 


Complex Bending 


The phenomenon of complex bending, which occurs 
when the resultant applied moment vector is not 
parallel to one of the principal axes, is quite common 
in aircraft practice. The term “‘complex’’ bending has 
been used instead of ‘‘unsymmetric’’ bending, since, 
even though the word unsymmetric is usually asso- 
ciated with the shape of the section, it is possible for 
the resultant moment vector to be nonparallel to a 
principal axis of a simple symmetric shape. The 
classic beam theory gives a formula for complex bending 
stress as follows:® 


a (myl, + mzI,,)x — (ml, + mylIry)y 
II, 2% (Izy)? 





fo (13) 

X and Y represent two mutually perpendicular 
centroidal axes; x and y represent the coordinate dis- 
tances from the center of gravity measured positive to 
the right and upward, respectively; m, equals the ap- 
plied moment about the X axis, the moment being 
positive when it tends to put compression in the 
“upper” fibers; m, equals the applied moment about 
the Y axis, the moment being positive when it tends 
to put compression in the “‘left’’ fibers, and J,, J, and 
I,, are the respective moments and product moment 
of inertia about the X and Y axes. 

This formula is applicable only in the elastic range; 
therefore, it needs to be modified for plastic bending. 
However, since the development of such a formula is 
beyond the scope of this article an alternative plan 
is proposed. In order to clarify the relationship be- 
tween section properties, applied moments, location of 
the neutral axis, and principal axes, a review of complex 
bending in the elastic range will be necessary for com- 
plete understanding of the alternative plan. 

Consider a rectangular section as shown in Fig. 8. 
The respective applied moment vectors about the 
centroidal X and Y axes, which in this case are principal 
axes, combine to give a resultant moment vector m at 
an angle of 6, with the X axis. Moment vectors are 
designated by double-headed arrows. They are to be 
interpreted by the left-hand rule (i.e., point thumb of 
left hand as indicated by the double-headed arrow and 
the natural curl of the remaining fingers of the left 
hand will designate the direction of moment). This 
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Fic. 8. 


implies, therefore, that the moment vector is the axis 
of the moment. 


Let a, = the angle the neutral axis makes with the 
X axis in the elastic range 


m,l. + mlz, (14) 


= tan 
ml, + myl,, 


for the general case but, due to the 
symmetry of this case, it will reduce 
to: 

tan—! (m,I,/m-,],) (15) 


Let 8, = the angle the resultant applied moment 
vector makes with the X axis 


= tan—'(m,/m,) (16) 


If the ratio of applied moments m,/m, is assumed to 
be approximately 5 and the ratio of moments of inertia 
I,/I, is assumed to be '/;5 when b/h = 4, then 


a, = tan~'(5/16) = 17.4° 
6, = tan-'5 = 78.7° 


The section and relative moments shown in Fig. 8 
represent the conditions of an actual lug that came up 
in practice. The largé difference in relative position 
of neutral axis and moment axis as indicated by the 
angle 8, — a, seems unusual. However, it has pur- 
posely been included in this discussion so as to familiar- 
ize the structural engineer with this phenomenon. The 
physical significance of this case indicates that, even 
though m, is small relative to m,, it is still capable of 
developing a relatively high compressive stress at 
point A on the section. In fact, the compressive stress 
in the upper fibers due to m, is not completely counter- 
acted by the tension stress in the upper right fibers 
due to m, until point B is reached. 

As long as the maximum stress/,, of the fiber fartherest 
away from the neutral axis (i.e., the upper left-hand 
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corner) is below the proportional limit of the material, 
the neutral axis location will remain as indicated 
above. 

On the other hand, as soon as f,, exceeds the propor- 
tional limit and the triangular distribution of stress 
relative to the neutral axis changes to a trapezoidal 
distribution, the neutral axis will rotate toward the 
resultant moment vector. This rotation will be a 
function of the relative distribution area on either side 
of the neutral axis as it rotates about the center of 
gravity of the section. Pending a further study of this 
phenomenon, it is proposed that the stress-ratio 
method® be applied with a straight-line interaction 
curve between the following stress ratios: 


R, = (m,/M,) and R, = (m,/M,) 


where M, equals the allowable moment about the 
principal axis XY and M, equals the allowable moment 
about the principal axis Y. 

The margin of safety may then be expressed as 


M.S. = [1/(R; + R,)] — 1 


Difference in Compressive and Tensile Properties of a 

Material 

All the preceding derivations and conclusions have 
been based on the assumption of identical tensile and 
compressive properties of the material. Although this 
may not always be true, the means of evaluating a 
new effective center of gravity of the cross section, 
based on a different distribution for the tensile and 
compressive bending stresses, are beyond the scope of 
this article. Since this situation also occurs in the 
present usage of the classic beam theory, it is proposed 
to neglect it, at least until the more important factors 
are taken care of. 


Combined Bending and Direct Load 


Many complications enter into a practical means of 
finding the position of the effective center of gravity. 
Even the simplest bending case in the plastic range will 
be subject to this effect when a direct axial load is 
included. Further study is definitely required with 
respect to this phenomenon. 

Here again it is recommended that the stress-ratio 
method® be used with an interaction between the three 
stress ratios represented by a plane surface: 


R, + R, + R, = 1.0 


where R, equals m,/M,, R, equals m,/M,, and R, 
equals f,/F,, and where f, is the applied axial stress, 
F, is the allowable axial stress, and x and y subscripts 
indicate ‘“‘about the principal X and Y axes,” re- 
spectively. 

The margin of safety in this case may be expressed 
as , 


M.S. = [1/(R, + R, + R,)] — 1 
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SUMMARY 


The following summary constitutes a presentation 
of the proposed method suitable for direct engineering 
use in aircraft work. 


General 


The conventional beam theory (based on the as- 
sumption that a plane section before bending remains 
plane after bending) gives a linear distribution of strain 
and stress in the elastic range (i.e., up to the propor- 
tional limit). In the plastic range, however, although 
the strain distribution is assumed to remain linear, 
the stress distribution corresponds with the stress- 
strain relationship for the material. 

An approximation of this distribution has been 
obtained which enables the prediction of the effects 
of shape and material properties on bending in the 
plastic range. 

Graphs are included (Figs. 9 through 17), each of 
which is applicable to a particular material. These 
graphs are to be used in connection with the following 
procedures. 


Simple Bending 


(A) Determine the ‘‘section factor” k = 2Q,,/(I/c). 
(Ifk > 2.0usek = 2.0.) Q, is the static (first) moment 
about the neutral axis of the area between the neutral 
axis and the extreme fiber. J/c ‘is the conventional 
section modulus. 

(B) Bending stress at the extreme fiber due to any 
applied moment, m. (1) Determine mc/I, (2) Enter 
the graph with mc/I and obtain the maximum stress, 
fm, at the extreme fiber for the proper ‘‘section factor”’ 
k. 

(C) Bending stress at an intermediate fiber due to 
any applied moment, m. 

(1) Obtain the stress at the extreme fiber from 
item (B). 

(2) Enter the stress-strain curve with this stress 
and obtain ‘the maximum strain e,, at the extreme 
fiber. 

(3) At the proper proportional strain e, of the 
intermediate fiber, which is at a distance y from the 
neutral axis, read the corresponding stress, where 
_ = (y/ C)Em. 

(D) Margins of safety in bending. (The margins 
must be determined from the bending moments since 
stress is no longer proportional to moment.) 

(1) Enter the graph with the allowable ultimate 
tensile stress F,,, to obtain an mc/I value for the proper 
“section factor’’ k. 

(2) Multiply this value by J/c to obtain the allow- 
able ultimate moment M,,. 

(3) Repeat (1) using the allowable yield stress 
F,, to obtain the allowable yield moment M,. 
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AM 265 HTA MAGNESIUM SAND CASTING 


The Effects of Shape and Material Properties 
on Bending in the Plastic Range 


Figure |7 


ress 
# att * 

in pst, Stress - Strain 
Curve 


F. 


Fx 107? 


tp 
Stress - Me Curves 


Me in psi x 107? 


0.6 0.8 


Stress-R Curve 


O16 
Strain, e 


Based on Minimum Guaranteed rties 


(4) The margin of safety on an ultimate basis: 
M.S., = (M,/jym,) — 1 


ju is an additional (multiplying) ultimate factor of 
safety, if any, otherwise use 7, = 1. m, is the applied 
ultimate moment. 

(5) The margin of safety on a yield basis: 


MS., = (M,/;,mz) — 1 


jy is an additional (multiplying) yield factor of safety, 
if any, otherwise use j, = 1.0. my is the applied 
limit moment and is usually = m,/1.5. 

This method of determining bending allowables has 
been substantiated by tests on aluminum-alloy castings 
and forgings, and on magnesium forgings of relatively 
thick-walled sections only. 

(E) Shear flow at any fiber due to any applied 
moment, m, and the corresponding shear, S. The 
formula that follows is applicable only in the special 
case where both principal axes are symmetric axes. 
This method of obtaining shear flow is included in 
order to indicate how the classic formula SQ/J is modi- 
fied by the nonlinear stress distribution. 

(1) Determine mc/I. 

(2) Enter the graph with mc/I and obtain an R 
value for the proper “‘section factor’ k, where R = 
dfo/dfm, the rate of change of intercept stress fy with 
respect to the maximum stress fp. 

(3) Substitute in the following: 





; 1 +R (46 — i) 
SQ Q 


q (shear flow) = —| ———— 
ITL1+R(R —- 1) 


where S equals shear load and Q equals static (first) 
moment about the neutral axis of the area A between 
the fiber in question and the extreme fiber. 


Local Crippling 


(A) Assume critical portion of section is developing 
its crippling stress. 

(B) Enter stress-strain curve with this crippling 
stress and determine the strain of the critical portion, é,. 

(C) Determine the proportionate strain of the ex- 
treme fiber e,,, where e¢, = (c/y)éy. 

(D) Enter the graph with this strain e,, to obtain 
an mc/I value for the proper “‘section factor”’ k. 

(E) Multiply this value by J/c to obtain the allow- 
able moment, M, which will subject the critical fiber 
to its crippling stress. 

For round tubes use the bending modulus of rupture 
curves available in ANC-5 since they include local 
crippling effects. 


Complex Bending 


This condition occurs when the resultant applied 
moment vector is not parallel to a principal axis. 

Let X and Y represent two mutually perpendicular 
centroidal axes; let x and y represent the coordinate 





% 
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distances from the center of gravity, measured positive 
to the right and upward, respectively; and let X’ and 
Y’ represent the principal axes. 

Moment vectors are designated by double-headed 
arrows and are to be interpreted by the left-hand rule 
(i.e., point left thumb in direction of vector and natural 
curl of fingers will designate the direction of moment). 


m, = the applied moment about the X axis, the 
moment being positive when it tends to 
put compression in the ‘‘upper’’ fibers 

the applied moment about the Y axis, the 


m, = 
moment being positive when it tends to 
put compression in the “‘left’’ fibers 

¢ = the angle the X’ principal axis makes with the 
original X axis, the angle being positive 
when measured in a_ counterclockwise 
direction. 

ye Y 

x 











Fic. 18. 


Any case of complex bending may be resolved into 
two cases of simple bending about the principal axes 
of the section. The principal axes are defined as 
mutually perpendicular centroidal axes about which 
the moments of inertia are a maximum and minimum, 
respectively, and about which the product moment of 
inertia is zero. 
oto &~4 

9 2 


z ‘ 
; a — sec 29 


where tan 29 = 2J,,/(I, — I). 
(A) Determine allowable moment M,, as described 
under simple bending. . 


Ome = 0.425+/1, A; (approximate formula) 
= LL? 
where A, = total area. “ee ee 


Cy = —xsing + ycos¢ 


where x and y are the coordinates relative to the original 
X and Y axes of the fiber fartherest away from the prin- 
cipal X’ axis. . 


ky = 2Q mz! /(I2+/C2") = < om” B = 
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(B) Determine allowable moment M, as described 
under simple bending, where 


Qmy 


cy 


0.425+/1,A , (approximate formula) 
xcosg+ysing 


Il 


where x and y in this case are the coordinates relative 
to the original X and Y axes of the fiber fartherest away 
from the principal Y’ axis. 


ky = 20 my/ (Ly /¢y’) 


(C) Determine applied moments about the principal 
axes. 


mM; = M,cos¢ + m, sin ¢g 
my = —m, sin g + m, cos ¢ 


(D) Determine stress ratios. 
Ry = my/M, and Ry = my/My 
(E) Determine margin of safety. 
M.S. = {1/[j(R: + Ry)]} - 1 


where j is the proper additional factor of safety de- 
pending on whether a limit or ultimate margin of 
safety is involved. 


Combined Bending and Axial Load 


(A) In a manner similar to that described under 
“Complex Bending’’ determine the stress ratios. 


Ry = my/My and Ry = my/My 


(B) Determine R, = f,/F,, where f, is the applied 
direct stress and F, is the allowable direct stress. 


(C) Determine margin of safety. 
MAS. = {1/[j(Rr + Ry + RJ} - 1 


where j is the proper additional factor of safety depend- 
ing on whether a limit or ultimate margin of safety is 
involved. 

(D) Secondary bending should be included when 
beam columns are being investigated. In such cases 
it is necessary to determine only the ultimate margin 
of safety. 
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Tabular Form 77 













































































































































































TABULAR FORM 77 - PLASTIC BENDING FACTORS BASED ON ss MATERIAL PROPERTIES’ AND SHAPE OF SECTION DATE Ps. 
REF. CODE: Oy * (Fepleed* + 3*O7@,+3=_ (where eg, strain at proportional limit) Ro. 
0 ® 9 KO) [a] GS] 
‘_ Val. from Stress-| First tercept] Values from Plot @ vs @| Plot @ versus these values 
tp” _ [Strom Curve | Mean | Mean [Interval] Moment] Geom. (Piet @ plot of @vs@| 4fm fo = 
Fey — 2 Strain | Stress Stress| Strain | Area |About NA) Factor |versus eo ~ Prnevened a eee Ue-ds ® + @-] 
y re 
Fmt] 8 Begin om we. 05[@n 05 [®, On: Qn + 1|©®@ ~~ T@rrer | Baos @ os @ here k equals: 
. EJP af mat |*@avi] *Onei] ~On] Tog] @ |-2@] * © |-0.| -@| © |-@ Fepw ye Ls] 20 
a. -. ey I 0.000 T | T 
= 2 a T mecal 
8 i 4, (Test) eel aa 
, ai 4 ae ae 
2 r Fy (MG) — : 
3 —- : 
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APPENDIX—ANALYTIC DERIVATIONS OF FORMULAS 


c 
m= of ofy dy 
0 


Moment of Resistance 


where z is width of section at point y (see Fig. 1 superimposed 


on Tabular Form 77). But 


f = fot Gm — fo)(y/c) 


2 2 
mm = 2 (for i - 4) dy 
A c c 


= fn(I/c) + fol2Qm — (I/c)] 
or 
mc/I = fm + folk — 1) (1) 
where k = 20,,/(I/c). 
Resisting Load P over a Portion of the Section 


P = load on that portion of the cross section between the 
extreme fiber and the fiber that is at a distance y from 
the neutral axis 


c 
P= fue 
y 


C 
f for dy + (In — fo 
y 


c 


( pa ) 
= pa + P—%o - 7.2 4 (a - 2 @) 
c c c 
A = area between the extreme fiber and the fiber that is at a 
distance y from the neutral axis 
Q = static (first) moment of area A about the neutral axis 


Shear Flow Distribution 


me 
— =fn + folk — 1) 
I 
cdm dfm . dfo 
-~—— = —— + —(k — ] 
Idx dx + dx. ) 
df m dfo Sc 
= —|1 —(k-—1)]/=— : 
dx | sa a.‘ | I (3) 
since S = dm/dx, from the relation P = fn(Q/c) + folA — (Q/c)): 
dP dfn Q , do dfm 4 
ao ae df. dx A me ie q (shear flow) 


Let R = rate of change of trapezoidal intercept stress fy with 
respect to maximum stress = dfo/dfm. 


_ Yn| Q _< 
¢™ dx E ? (A °)] (4) 


But from Eq. (3): 
dfm/dx = (Sc/I)/(1 + R(k — 1)] 
Substitute in Eq. (4): 


_ (Sc/D{(Q/c) + RIA — (Q/o)}} 
=” 1 + Re — 1) 


Therefore, 





Factor out Q/c 


_ SQf1+ R[(Ac/Q) — 1] 
:= t 1 + R(k — 1) 





(5) 


Sesenations Factor to Reduce Test Data to Minimum Guaranteed 
‘alues 


In Fig. 2 (superimposed on Tabular Form 77): 
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f = Fu + (Fu — Fry)/cly 
a Fy’ + ((Feu’ — Fry’) /cly 


Fy, — F, 
fotes fo Pe + yim Flay 
Correction _ Jo 2 re ¢ we 
factor a 
7] Fin’ =e F, , 
Sx ” >| Fe’ + ya = Fay 
0 


c a 
Fu + (Fru — Fouls 








c mr 
Fy’ 5 + (Fis! — Fry’); 
(Fry + 2Fru)/(Fry’ + 2F ru’) (6) 


Approximation to Qm 
Let 20m = KpA:. K is a relationship factor, p is radius of 


gyration, and 4A; is total area. 
Type (a) rectangular cross section (Fig. 3 superimposed on 


Tabular Form 77): 


= y.. - Pe moe? 
At bh 3.46 
pA, = bh?/3.46 
20m = (2bh/2)(h/4) = bh?/4 
Therefore, 
K = 2Qm/A:t = 0.865 


Type (b) diamond cross section (Fig. 4 superimposed on 
Tabular Form‘ 77): 


exe — = 4/245 


(2bh/2) 
pA, = bh®/2.45 
20m = (2bh/2)(h/3) = th?/3 
_ (bh®/3) 
K « ——__-_ = 0317 
(bh?/2.45) wel 


Type (c) hour-glass cross section (Fig. 5 superimposed on Tabu- 
lar Form 77): 





_ . [(2bn2/4) SeuEre 
| = qe V h?/2 = h/1.414 
| pAy = bhY/1.414 
20m = (2bh/2)(2h/3) = bh?/1.5 
ao 


bh?/1.414) — 
A general value of K = ().85 has been assumed, which in most 
cases is conservative. 








Relative Values of Trapezoidal Stress Intercepts 


In Fig. 6 (superimposed on Tabular Form 77): 

Let Q, designate the actual static moment about the N.A. 
of the area OABCD. 

The equivalent trapezoid OABE has a trapezoidal stress inter- 
cept QE = fo and a maximum stress AB = f,, at the extreme fiber. 

The equivalence of the trapezoid lies in the fact that /p has 
been determined so that the static moment of the trapezoidal 
area about the neutral axis is equivalent to Q,. 

This may be written as 


Q; fo c#/2 + fas m5 fo) c?/3 (7) 
= (foct/6) + (fme?/3) — 


from which 





fo = (6Q./c*) — 2fm (8) 
(Column 8 on Tabular Form 77) 

The general expression for equivalent first moments about 
the neutral axis for a stress-strain distribution and a trapezoidal 
distribution should include the variation of the cross-sectional 
width z from the neutral axis to the extreme fiber. This may be 
written as 


c 
m = 2 {yeh ay (9) 
0 


° 
When z is a constant as in the rectangular type of cross section, 
the above reduces to 


m = 22 {x ay = 220, (10) 

/0 
which indicates that the resisting moment m is directly propor- 
tional to the stress-strain static moment Q, and, therefore, to 
its equivalent trapezoid as represented by OABE in Fig. 6 super- 


imposed on Tabular Form 77. 
However, in the case of the diamond-type cross section where 


b is the width at the neutral axis: 
z= b — y(b/c) (11) 


c 
m=2 [ vefay 7 2 fy — y(b/c))f dy 
0 


2 . 


260; — (6/c)Is (12) 


since J,, the moment of inertia about the neutral axis for area 
OABCD and a symmetric area below the axis, is 


iC 
I, = 2 {yay (13) 
0 


oe 


This expression for moment may be written in terms of a true 
equivalent trapezoidal stress intercept designated as foa and the 
maximum fiber stress f,,: 


I 
m = —[fm + foa(k — 1)] (see Eq. 9) 
c 


Since k = 2 for a diamond section (see Table 1): 


m = (I/c)[fm + foa] (14) 
Substitute Z/c for a diamond section = bc?/6 and equate to 
Eq. (12): 
bc? b 
m = —(fm+ foa) = 200, — -I, (15) 
6 c 
120. 61, 
fea ae Se at my (16) 
c c 


But QO, = fo(c?/6) + fm(c?/3) (exact value) and J, = (foc?/6) + 
(fmc?/2) (approximate value). Therefore, 


Bi Af «: “) Cf sf « “) 
2 \ iq ting ca og + Img] fn 


fo (approximately within 2 per cent) (17) 


Soa 


In the case of an hour-glass section where 5 is the width at the 


extreme fiber: 
z = (b/c)y (18) 


€ c 
m= 2f wfdy =2f y0/oylfay = W/o a9) 
o 


This expression for moment may be written in terms of a true 
equivalent trapezoidal stress intercept designated as fy, and the 
maximum fiber stress fm; 


m = : [fm + forn(k — 1)] (see Eq. 9) (20) 











CO SSS eee 

















BENDING STRENGTH IN THE PLASTIC RANGE 151 














































































































i k = */;f hour-glass section: “a fuc® 
vince fo San ein Deengane eet Substitute 7, = ms +° = (approximate value). 
T ) < 
m= 1 fn +] 6 3 3 
fon = — E ; + fm 4 — 3fm 
Substitute J/c for an hour-glass section = bc?/2 and equate ; ' i 
to Eq. (19). = fo (approximately within 2 per cent) (23) 
bc? fon # Since true values of the trapezoidal stress intercepts for the 
welder” Sm + me ag (b/c). (21) three basic types of sections are within 2 per cent of each other, 
, the resisting moments computed on the basis of the trapezoidal 
ha te 6 — (22 stress intercept fo for a rectangular section are considered satis- 
° eo mn 5 factory 
a.  TaBLe 4 : : . 
Test Data‘ 
1|2 3 4{|5/ 6 | 7 | 8 9 10 11 | 12 13 14 15 16 17 
ai. mo Predicted! Actual [Correo /Actual/Actual |, 2atio laverage 
g Drewing | * | 2 |section Properties I Allowable} Physical tion Moment | Moment |», pre. X-Ray 
= and ql I k lat Fy, |Ultimate | Properties | Factor fram |Correct-gicted |Rating 
$13 Part & s > 2 lnax (7) for Moment From Test Test ed Moment | for 
~ | ee) tu* Fey ing 2° 
X—5495-2 | S31 | 4/0.1800/0.2816)1.564} 95,700} 17,200 | 68,955/60,050) 0,910 20,860}18,970 | 1.104 
=-2/ 80} 4/0.1800/0.2816/1.564} 95,700] 17,200 | 68,955/60,050} 0.910 20,900}19,000 | 1.105 - 
-3 | 33 | 4/0.1314/0.1793/1.365| 86,600; 11,379 |68,955/60,050) 0.910 13,480/12,250 | 1.079 - 
g -4/13 | 2/0.1239/0.1642|/1.327| 84,500/ 10,470 | 66,700/56,950) 0.945 11,800/11,130 | 1.064 - 
89 4/14] 2/0.1239/0,1642/1.327| 84,500) 10,470 | 66,700/56,950) 0,945 11,840/11,200 | 1.076 - 
3 -5| 35 | 3/0.1165/0.1494/1.283/] 82,000} 9,553 | 68,650/59,450) 0,915 10,980}10,000 | 1.045 - 
as -6 | 37 | 3/0.0971/0.1296/1.335| 85,000} 8,254 | 68,650|/59,450| 0,915 9,360) 8,560 | 1.039 ~ 
Q 3 X=309S-1 | 1] 1/0.1667/0.2500/1.500/ 95, 15,800 |67,940|59,975| 0.920 18,120) 16,670 | 1.054 om 
ola -1;/ 2] 1/0.1667|0.2500/1.500} 95,000] 15,800 |67,940/59,975|) 0.920 17,680/ 16,250 | 1.c28 - 
e -2/ 3] 1/0.1603/0.2240/1.397| 88,800] 14,235 | 67,940/59,975) 0.920 15,760/ 14,500 | 1.020 ~ 
B -2| 4] 1/0.1603/0.2240)/1.397/| 88,800; 14,235 |67,940/59,975| 0,920 16,400/15,070 | 1,059 ~ 
7 -3/| 5] 1/0.1260/0.1600/1.270| 81,200; 10,231 /|67,940/59,975) 0.920 11,440} 10,530 | 1.030 - 
wid -3| 6] 1/0.1260/0.1600/1.270; 81,200} 10,2351 | 67,940/59,975) 0.920 11,560} 10,650 | 1.040 - 
-1/ 7] 2/0.1667/0.2500/1.500| 56,250) 9,377 |50,715/35,125) 0.792 12,000} 9,500 | 1.013 - 
§ -1;) 8] 2/0,.1667/0.2500/1.500| 56,250} 9,377 |50,715/35,125) 0.792 12,320) 9,750 | 1.040 ~ 
ae -2/ 9] 1,0.1603/0.2240/1.397| 53,100| 8,512 |50,725|32,790| 0.805 11,160} 9,000 | 1.057 - 
a =2| 10] 2/0»1603/0.2240/1.397| 53,100} 8,512 /|50,715|35,125) 0.792 11,420) 9,050 | 1.062 - 
2 -3| 11] 1/0.1260/0.1600/1.270|) 49,400} 6,224 |50,725/32,790; 0.805 8,120] 6,540 ; 1.050 - 
4 -3|12/ 1/0.1260/0.1600/ 1.270} 49,400} 6,224 |50,725|/32,790) 0.3805 7,840] 6,310 | 1.013 - 
a ¥ -1/ 48 3/0.1667/0.2500/ 1.500} 44,000; 7,335 |35,9350/27,500| No 8,160] 8,160 | 1.114) Fair 
ee -2| 49] 3/0.1603/ 0.2240) 1.397) 41,500} 6,652 [35,930/27,500| Reduction| 7,552) 7,552 | 1.135/ Poor 
33 -3 | 75 | 3/0.1260/0.1600/1.270| 38,750} 4,883 |37,060/28,675| Factor 5,976] 5,976 | 1.225 7 
nm -3 | 78 | 3/0.1260/0.1600/1.270} 38,750; 4,883 /|37,060/28,675 5,440) 5,440 | 1.113 - 
8] | X-5252-1/ 16 | 2 | 0.0982) 0.1670) 1.700} 107,000] 10,507 |66,700)56,950) 0.945 12,480| 11,800 | 1.123 - 
Sle D -1/ 79} 4 | 0.0982/0.1670/1.700| 107,000] 10,507 |68,955/60,050} 0.$10 12,280} 11,170 1.063 ° 
° <8 -2] 17] 2 | 0.0854/0.1310/1.535} 97,100; 8,292 |66,700|/56,950| 0.945 9,880] 9,340 1.125 - 
a| -3/ 19] & | 0.0580/0.0810}1.396] 88,500} 5,133 /|68,975/59,700) 0.915 5,660} 5,150 1.002 
3 © & -1| 39] 3 | 0.0982/0.1670/1.700] 49,000| 4,812 |35,930/27,500| No 6,424] 6,424| 1.336 | Feir 
alae -2| 40] 4 | 0.0854/ 0.1310)i.535] 44,750| 3,832 |31,935]26,900|Reduction| 4,328] 4,328] 1.131 | Feir- 
ae -3] 41] 3 | 0.0580/0.0810/1.396] 41,400] 2,401 |35,930/27,500| Factor | 3,456) 3,456 | 1.440 |Feir 
o 
wo | X-5496-2] 29} 2 | 0.1176/0.1806/1.537| $7,200] 11,431 |66,700/56,950/ 0.945 (|13,040/ 12,320] 1.078 - 
sled -3] 27] 5 | 0.0950|0.1378|/1.453| 92,200| 8,759 |68,975|59,700| 0.910 9,980] 9,080] 1.037 * 
g/S5 -4| 21] 5 | 0.0884/0.1212/1.372| 87,300| 7,717 |68,975|/59,700| 0.910 8,200] 7,450] 0.965 . 
. f -5| 25] 5 | 0.0685/0.1020/1.490] 94,400} 6,466 |68,975/59,700| 0.910 7,780} 7,080 1.0¢5 ~ 
& -6| 23] 3 | 0.0605/0.0836|1.383} 88,000] 5,324 |68,650/59,450| 0.915 5,960] 5,450| 1.024 . 
et X-3494-2]) 62] 2 | 0.1296/0.1838/1.419/ 42,000] 5,443 |36,295/28,725 5,740} 5,740 1.054 | Poor- 
g ee -3] 60] 2 | 0.1082/0.1444/1.337] 40,400] 4,371 |36,295/28,725| No 4,720| 4,720] 1.080 |Poor+ 
s sz -4/ 61] 2 | 0.0997/0.1259/1.264) 38,600 3,848 36 ,295| 28,725 |Reduction| 4,752} 4,752 1.235 |Fair- 
S/o 2 -5| 59| 2 | 0.0711/0.0975/1.372| 41,000| 2,915 |36,295/28,725| Factor | 3,136| 3,136 | 1,075 |Poor+ 
o -6|63-A 4 | 0.0660/ 0.0864/ 1.310; 39,600 2,614 31,935] 26 ,900 2,920 2,920 1.117 |Poor 



























































A New Method of Longitudinal Control for 
Aircraft by Use of an Adjustable Angle of 
Attack Balance 


ROLAND J. WHITE* 
Curtiss-Wright Corporation 


INTRODUCTION 


Cn PRESENT-DAY TYPE of large airplane has ap- 
proached a definite limit in two respects: first, it 
must be loaded over a restricted center of gravity range 
(c.g. range); second, as the airplane size and gross 
weight are increased, the control forces become too 
great unless properly reduced, which so far has been a 
troublesome problem. 

The c.g. range of the airplane has the forward limit 
fixed by the maximum control force that can be pro- 
duced by the horizontal tail to meet a given landing or 
maneuvering requirement and the rearward c.g. limit 
established by the ability of the airplane to meet cer- 
tain stability requirements. Lately, these stability 
requirements have offered much discussion, mainly be- 
cause airplanes may still be flown at c.g. positions 
further aft than the position generally limited by the 
requirements. 

In this paper, the author presents a new method of 
aircraft longitudinal control which it is believed will 
permit the c.g. range of the airplane to be substantially 
increased, as well as reduce the control forces in a more 
satisfactory manner than is customary. This method 
essentially provides a cockpit control by which the pilot 
may regulate the elevator hinge moment characteristics 
according to an estimate of the airplane c.g. and enables 
the pilot to make the airplane stable during a climb or 
in level cruising flight or maneuverable when he so de- 
sires. The cockpit control is a means for regulating 
the amount the elevator hinge moment changes for a 
given change in angle of attack of the airplane horizon- 
tal tail as expressed by the elevator hinge moment co- 
efficient derivative 0Cy,/Oa; = Creag: 

To obtain this regulation a method is proposed where- 
by the forward part of the stabilizer is made movable 
and is connected to the elevator by an adjustable link- 
age system such that it may move in a regulated man- 
ner with the elevator. A second method consists of 
allowing a small surface within the stabilizer to move 
in an adjustable relation to a tab attached to the ele- 
vator, the internal surface being vented to the airfoil 
surface. 
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In the first part of the paper the elevator control 
forces are analytically calculated in order to show how 
the various physical characteristics of the airplane effect 
the longitudinal control forces. Following this, an ex- 
ample airplane having arbitrary dimensions is selected, 
and control forces are calculated in order to demonstrate 
the use of the developed equations and to indicate the 
improvements that may be realized if Cy,, is properly 
regulated. In the second part of the paper wind-tunnel 
results of the two methods of regulating Cy,, are given 
along with charts and equations permitting the engineer 
to design properly such a control system. 

It is believed that future airplanes may have their 
c.g. range extended to twice the present permissible 
range if proper regulation of the hinge moment deriva- 
tive Cyeq is made possible. 


SyMBOLS Useb 


b = wing span (ft.) 

ly = distance from c.g. to elevator hinge line (ft.) 

MAC = wing mean aerodynamic chord (ft.) 

h = length of control column (ft.) 

o = angular position of control column 

ce = effective moment arm of internal pressure sur- 
face (ft.) 

Ces 08 = MAC of elevator and tab, respectively (ft.) 


D,, D2, Ds = airplane aerodynamic constants defined by Eqs. 
(6), (10), and (13) 


Cc = constant involving the mechanical advantage of 
the elevator control system 
R = ratio of moment transmitted from balance sur- 


face to tab or elevator to that developed by the 
balance surface (positive if moment tends to 
deflect elevator or tab in a positive direction) 


Ry = wing aspect ratio 

Ry, Rn = aerodynamic constants defined by Eqs. (25) and 
(28) 

N = speed ratio given in Eq. (21) 

o = density ratio 

V = air speed (ft. per sec.) 


V = tail volume = (S,/S)(1,/ MAC) 

tail efficiency factor = (V at tail/V)? 

angular pitch velocity (rad. per sec.) 

airplane load factor 

= elevator control forces given by Eqs. (23) and 

(27) 

G = distance from airplane aerodynamic center (a.c.) 
for tail-off to airplane c.g. in fraction of wing 
MAC (positive tor c.g. forward of a.c.) 

| = airplane wing loading (lbs. per sq.ft.) 

S, S; S., = areas of wing, tail, elevator, internal balance 

So, Ss surface, and tab, respectively (sq.ft.) 
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Cr = airplane lift coefficient (neglecting tail lift) 

Cro = airplane lift coefficient such that C,, tail-off equals 
C,, tail-on for e = 0 

Crm = maximum lift coefficient 

Cr = tail lift coefficient 

k., ke = elevator and external balance effectiveness 
(taken as positive) 

Map, Mt = dC,/da for wing and tail (C,/deg.) 

ee = airplane pitching moment coefficient 

ts = airplane pitching moment coefficient at C, = 0 
fore = 0 

Gs = airplane pitching moment coefficient due to the 
horizontal tail 

Cnat. Cna = airplane pitching moment coefficient developed 


due to angular pitch velocity for tail and for 
airplane minus tail 

Cmg = O0Cm/0O(qb/2V) for airplane minus tail 

Cue, Cus = hinge moment coefficient of elevator and tab 
based on their respective chords and areas 
(a line over the ‘‘C’”’ indicates the elevator or 
tab acting with balance surface) 


Cun = hinge moment coefftcient of external balance sur- 
tace based upon the elevator chord and area 
CHa = OCyn/OCzt for e = 0 


Cype, Cupe= OCyn/O8, OCyn/Oe for Crt = 0 
Cree CHea,= OCye/Oe, OCye/Oa, for a or e constant. The 
etc. second subscript indicates the variable of the 
derivative 
P = ratio of pressure difference between the upper 
and lower surfaces of the airfoil at some chord- 
wise position (x/c) divided by (p/2)V? 


Basic CONTROL Force EQuATIONS 


In developing the basic control force equations the 
following assumptions will be made: (1) Airplane pitch- 
ing moments and their components are linear with C,; 
(2) Horizontal tail lift and elevator hinge moments are 
linear with tail angle of attack and elevator deflection; 
(3) All accelerations are neglected except the normal 
acceleration due to a curved flight path. 

Assumptions (1) and (2) are convenient but may not 
be applicable to many present airplanes; however, it is 
believed that for future large airplanes these conditions 
must be realized before efficient longitudinal control will 
be obtained. 

Referring to the list of symbols, the airplane pitching 
moment coefficient in recilinear flight for the airplane 
without tail may be written as 


Cai = <n = GC, (1) 


To this pitching moment coefficient must be added 
that due to the tail. The tail pitching moment coeffi- 
cient is influenced by the down-wash from the wing, 
which may be represented by da,/da. From elemen- 
tary down-wash theory for the tail a normal distance 
behind the wing, 


da,/da = [1 — 2(57.3)m,/rR,] (2) 
The tail lift coefficient is expressed as 
Cr: = ma, + kee) (3) 


and defining the tail volume as 
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V = (S,/S)(l,/ MAC) (4) 
the tail pitching moment coefficient may be written as 
Cnt = — mVD(Cr — Cro) — 1Vmke (5) 

where 
D, = (m,/m,)da,/da (6) 


Having determined the airplane pitching moment 
components for rectlinear flight the effects of a curved 
flight path will be considered. For a pull-up from level 
flight the airplane lift coefficient is 


C, = nl,/(p/2) V? (7) 
and for a given load factor m the angular pitch velocity 
q, expressed in radians per second, is 

q = (nm — 1)g/V (8) 
The above relations permit the increment of angle of 
attack at the tail due to g to be calculated as 

Aa, = [(n — 1)/n|D2C, (9) 
where 
Dz = 2.19l,0/+/ 01 le (10) 


here, 2.19 = 0.002378(32.2)(57.3)/2. The correspond- 
ing damping moment in pitch because of the horizontal 
tail then becomes 


Cosas = niVmDel(n = 1) ‘n\Cr (11) 


In addition to this damping moment the wing and fuse- 
lage will contribute a damping moment that is about 25 
per cent of that of the tail for a normal type airplane. 

In the present notation this damping moment coeffi- 
cient may be expressed as 


o—_ = i(u - 1) n\D3Cr (12) 
here 
Ds = 0.0191 Cngob/Lip (13) 


where 0.0191 = 0.002378(32.2)/4. Now combining 
the ‘pitching moment increments from Eqs. (1), (5), 
(11), and (12) the total airplane pitching moment equa- 
tion becomes, 


Cu = [Cro + niVD:Cz0) a |G + n.VD, + 
(=—+) (n.VmDe fit Dd») |e, _ n,Vm.ke (14) 
n 


In order to find the most forward value of G for trim of 
the airplane let Ch, = Cpm, Cn = 0, Re = (Ree) m; then 


-a]-[n(-€) 
Ga ( — ~ A1 ~ of DB [1 — 2) + 
J le ’ : ; Com 


(n — 1) 


(m (Re) m/Crm) + ~ Dam. | (15) 





This forward c.g. position represents the limit for con- 
trol, assuming, of course, that the pilot can exert 
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enough control force to produce the maximum elevator 
deflection desired. The next problem is to relate the 
elevator control force with the airplane pitching mo- 
ment coefficient. To do this, express the elevator hinge 
moments by the equation 


Cre = Cread: + Creel 


During a maneuver a definite value of e will be required 
which may be found from Eq. (14) for the condition 
that C,, = 0; also a corresponding angle of attack of 
the tail will occur which may be determined from Eqs. 
(3) to (9). It will now be essumed that the airplane 
is in trim in level flight at a speed V; giving Cy, and a 
maneuver is made to a new value of C, = Cy2; then the 
change in elevator angle necessary to produce this 
change in C, as found from Eq. (14) is 


(16) 


—1/ G D 
Ae = “Us + Di + mD, — >| (Crz — Cr) 
m ken V nib 
(17) 
and the corresponding change in tail angle of attack is 
[(D,/m,) + D2)(Crz — Cr) (18) 


From these relations the change in elevator hinge 
moment coefficient becomes 


ACiea baad Cue: Aa, —- Cree Ae 


Aa, = 


(19) 


The elevator control force F will depend upon the me- 
chanical advantage of the control system which is as- 
sumed constant for all control positions and indicated 
by the derivative d¢/de. As the elevator hinge mo- 
ments are based upon the elevator area and chord, the 
maximum control force required to maneuver is 


F= [.S.c./h(do/de) | [(p/2) V?)( ACu.) 


There are two types of manuevers to be considered. 
The first type assumes that C;2 is obtained by allowing 
the air speed to change to V2 with no appreciable accel- 
For this case let nm = 1 and introduce 


N = Vi/ Ve 


(20) 


eration. 
(21) 
then 

Crz — Cr = [lw/(e/2)Vi7][N? — 1] (22) 


Using the relations of Eqs. (17), (18), (19), (20), and 
(22) this control force, denoted as Fy, may be expressed 


as 





Fy, = CRy[1 — 1/N?] (23) 
where 
C = Sdy/h(do/de) (24) 
| ( Cae G | x 
Ry = — D,{1—k, —*) + -—> 3) 
: M (Re ; Cree mV { >) 


The second type of maneuver is produced by a pull-up 
or highly banked turn that produces an acceleration but 
does so without an appreciable change in air speed. 
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For this case the velocity V of Eq. (2) remains constant 
at V = V, and 


Cre — Cu = (n — 1)l./(p/2) Vi? (26) 


Using this relation in a similar manner, as was done in 
the first case the control force, denoted as F,, that is re- 
quired to perform this maneuver will be independent of 
the air speed and given by the equation 








F, = CR,(n — 1) (27) 
where 
RK, = 
—Gnel (p, + Dyn, (1 —k Gx) 4 (2?) ] 
M (Re Hee mV J 


(28) 


Eqs. (23) and (27) represent the basic control force 
equations and may be used in combination if necessary 
to calculate any longitudinal control maneuver. 

The interpretation of these basic control force equa- 
tions may best be studied by applying them directly to 
an airplane. This will be done in the following section. 


CALCULATIONS AND DISCUSSION OF CONTROL FORCES 
FOR A TYPICAL LARGE AIRPLANE 


A typical large airplane of arbitrary dimensions has 
been chosen and the control forces have been calculated. 
It will be assumed that the elevator has a high degree 
of aerodynamic balance (Cy,. = —0.002) and can main- 
tain this balance at all elevator angles from e = +10° 
toe = —25°. A constant mechanical advantage of the 
control system will be assumed such as to permit an 
18-in. fore-and-aft movement of a 3.0-ft. control col- 
umn. The horizontal tail is assumed to have a 40 per 
cent chord elevator that is capable of producing a maxi- 
mum value of (k,€)m = —10 which is used in deter- 
mining the forward c.g. position for trim. This as- 
sumes that k, will reduce in value at a large elevator 
angle; hence, the control force required to produce 
e = 25° elevator deflection will be larger than those 
calculated. This condition would require a special in- 
vestigation and will not be made here. 


Airplane Dimensions 


Gross weight = W = 50,000 lbs. 

Wing area = S = 1300 sq.ft. 

Wing loading = /,, = 38.5 lbs. per sq.ft. 
Wing span = b = 102 ft. 

Wing MAC = 13.35 ft. 

Wing aspect ratio = R = 8.0 

Slope of wing lift curve = m, = 0.080C,/deg. 
Distance from c.g. to elevator hinge axis = 1, = 40 ft. 
Tail area = S, = 175 sqft. 

Tail area ratio = S,/S = 13.45 per cent 
Overall span of tail = 28.0 ft. 

Tail aspect ratio = 4.48 
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Elevator chord = 2.50 ft. 

Elevator area = 65 ft. 

Slope of tail lift curve = m, = 0.050C,/deg. 

Elevator hinge axis 60 per cent chord position from 
L.E. of stabilizer 

Elevator effectiveness = k, = 0.670 

Maximum effective change in tail angle of attack due 
to deflecting the elevator = (ke), = —10° 

Tail efficiency = (V,/V)? = 0.90 

Effective tail volume = 0.90(0.1345)(40/13.35) 
0.363 


Calculated Constants 


Cn = —0.006 

Crm = 1.5 

a.c. = 0.15 MAC from L.E. for tail-off 

D, = 0.397 D; = 0.010 (Estimated) 


D, = 2.40 for ¢ = 1.0 C = 2540 


These dimensions are first used with Eq. (23) to cal- 
culate the control force experienced during an unacceler- 
ated pull-up, assuming that the airplane is in trim at a 
given air speed and that a climb is started such to reduce 
the air speed by 20 per cent before a change in trim is 
made. These forces are calculated for various c.g. 
positions using different values of Cy,. and plotted in 
Fig. 1. From Eq. (15) the most forward c.g. for trim is 
15.5 per cent of the wing MAC from its leading edge. 
It will now be assumed that in obtaining an aero- 
dynamic balance giving Cye-e = —0.002 a value of 
Cuea = —0.0005 will result for the normal airplane. 
If, however, Cy., is regulated such that dCm/dC, for 


o——-—— NORMAL A/RPLANE Creu” 0005 
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Fic. 1. Calculated control forces for example airplane. 
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the airplane with elevator free maintains a constant 
value of static stability—-say, —0.05 for the elevator 
free condition—then the curve labeled 5 per cent static 
stability in Fig. 1 will result. If Cy., could be con- 
trolled to +0.0020 for this airplane a c.g. range from 
15.5 to 34.5 per cent of the MAC would result while 
still having 5 per cent static stability. In this case 
the most rearward c.g. for dCm/dC, = 0 for elevator 
free would be 39 per cent of the WAC. IfdCm/dC, = 0 
for elevator free is adopted for the criterion in deter- 
mining the most rearward c.g. position, the normal air- 
plane would have a 11 per cent c.g. travel, while with 
proper control of Cy,, a c.g. range of 23.5 per cent would 
result. It is well to note here that had the value of 
Cee = —0.002 been obtained by using a normal type 
balancing tab a value of Cy,, = —0.004 would result 
which would make the airplane unstable at all c.g. posi- 
tions. 

In Fig. 1, curves of the control forces experienced dur- 
ing a 3g pull-out at constant air speed are also shown as 
calculated from Eq. (27). During this type of ma- 
neuver the most forward c.g. position is limited to 19 
per cent if it is necessary to maneuver at a reduced 
speed such as to give C,,, = 1.5 at a 3g maneuver. 
Now if the most rearward c.g. is selected as in 
the above paragraph for the wunaccelerated flight 
maneuver, then the normal airplane would have a c.g. 
range from 19 to 26.5 per cent which would result in a 
maximum control force of 90 Ibs.; while, if Cy. could 
be controlled such as to maintain zero static stability 
for elevator free, the airplane would have a c.g. range 
from 19 to 39 per cent with a maximum control force of 
55 Ibs. This indicates that by regulating Cy,, the c.g. 
range has been extended 2.67 times while reducing the 
maximum control forces by 39 per cent. It is felt that 
if the airplane would have an adjustable angle of attack 
balance that would double the c.g. range it would be 
sufficiently satisfactory. 

Further interpretation of these figures will not be 
made, but they show the resulting advantage that may 
be obtained if proper regulation of Cy,. is made during 
flight. By proper use of the basic control force equa- 
tions the control forces for any proposed airplane may 
be calculated and the requirements for the detailed 
aerodynamic design established. 

There are other factors that may limit the c.g. range 
indicated by the above equations and the extent to 
which Cy, may be controlled. These factors cannot be 
discussed here but will be outlined as follows: 

(a) Effect of landing with and without flaps. 

(b) Effect of k, being reduced at large elevator angles 
including trim tab effect. 

(c) Effect of dynamic stability of the airplane for the 
elevator free condition. 

(d) Effect of slipstream on the airplane pitching mo- 
ments. 

(e) Effect of compressibility at high speed. 
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Fic. 2. Rear quarter view of external elevator balance test setup 
(University of Washington wind tunnel). 
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EXTERNAL ELEVATOR BALANCE 


INTERNAL TAB BALANCE 
Fic. 3. Diagrammatic cross-sectional view of angle of attack 


balances tested, omitting for simplicity the means of adjusting 
the linkage. 
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DESCRIPTION OF TESTS MADE ON Two TYPES OF 
ANGLE OF ATTACK BALANCES 


External Elevator Balance 


The first type of angle of attack balance was tested in 
the University of Washington wind tunnel in Seattle 
while the author was connected with this University. 
The model consisted of a N.A.C.A. 0009 airfoil having 
the forward portion hinged at the 30 per cent chord axis 
to form an external elevator balance surface. This 
surface was connected by mechanical means to a 35 per 
cent chord elevator such that a downward movement of 
the elevator caused a small upward angular movement 
of the external balance surface, which in effect meant 
that a fraction R of the hinge moment developed by the 
balance surface was transmitted to the elevator. A 
rear quarter view of the test setup showing the wire pul- 
ley arrangement for measuring elevator hinge moments 
is shown in Fig. 2. In Fig. 3 a diagrammatic sketch of 
the typical linkage is shown, but for simplicity the 
means of adjusting the links is not shown. The hinge 
gap of the balance surface was sealed at the center by 
a squeezed rubber tube, and a gap of 0.10 in. occurred 
at the airfoil surface for the balance surface in the un- 
deflected position. The elevator gap was not sealed but 
was made small. Model dimensions are given in 
Table 1. 
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Fic. 4. Measured lift and hinge moment coefficients with 
external elevator balance. 
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Fic. 5. Comparison of measured and calculated hinge moment 
coefficient derivatives for external elevator balance. 


Some results of these tests are shown in Figs. 4 and 5. 
In Fig. 4 the elevator hinge moment, coefficient Cue for 
a, = 0, is plotted against elevator angle for different 
values of R. These curves are linear up to e = 15° 
after which their slopes are constant. The lower 
curves show how the lift coefficient varies when the 
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TABLE 1 
Model Dimensions 


Internal 





- External 
Elev. Bal. Tab Bal. 
Model Model 

Wing span 36 in. 60 in. 
Wing chord 15 in. 24.54 in. 
Wing aspect ratio 2.4 2.79 
Diameter of circular end plates 24 in. 36 in. 
Elevator chord aft 5.25 in. 8.07 in. 
Tab chord 2.0 in. 
Tab span 42.0 in 
Balance effective moment arm or chord 4.50 in. 8.0 in 
Balance effective area 162 sq.in. 57 sq.in 





elevator is allowed to float downward. On this same 
curve the lift coefficient for e = 0 is plotted for compari- 
son, along with the curve for a normal free elevator 
(R = 0). In Fig. 5 values of the hinge moment deriva- 
tives for different values of R are shown plotted along 
with calculated values of these parameters determined 
by the method presented in a later section of this 
paper. The maximum lift coefficients are not shown, 
but it was found that the effect of the movable nose 
balance for the value of R tested did not reduce the 
maximum angle of stall to any appreciable extent. 


Internal Tab Balance 


The second type of angle of attack balance was 
tested at the Massachusetts Institute of Technology 
and consisted of a N.A.C.A. 0010 airfoil having a 40 
per cent chord elevator with a 20 per cent elevator 
chord tab having a span of 70 per cent of the elevator 
span. This tab was connected to an internal tab bal- 
ance surface located inside the stabilizer in a typical 
but somewhat different manner to that shown in Fig. 
3. The internal balance surface had a slack diaphragm 
seal with the regions on each side vented by a series 
of holes located in the metal airfoil surface. An ad- 
justable means was provided to change the movement 
of this internal surface in relation to the tab movement. 
Elevator hinge moments were measured by use of elec- 
tric strain gauges, and a fractional horsepower electric 
motor was provided to change the elevator angle with- 
out stopping the wind tunnel. Two autosyn indica- 
tors were used to measure the elevator and tab angles. 
Fig. 6 shows a detail view of the internal balance sur- 
face and autosyn angle transmitters. In this model the 
elevator, tab, and internal surface were carefully stati- 
cally balanced. Model dimensions are given in Table 
Z. 

Typical test results are shown in Figs. 7 and 8. In 
Fig. 7 the lift coefficient plotted against angle of at- 
tack for the elevator free is shown for the internal bal- 
ance surface vented at 10 and 15 per cent of the airfoil 
chord for R = —0.0426. This shows that the tab 


effectively controls the elevator up to the limit of the 
tab which was near 30°. 


The maximum tab deflection 





FOR AIRCRAFT 15 





Detailed view of internal tab balance showing elevator 
and tab angle autosyns. 
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Fic. 7. Elevator free tests of internal tab balance for pressure 
vents at 10 and 15 per cent of chord (R = —0.0426). 


varied to some extent because of restraints on the me- 
chanical mechanism which changed for different eleva- 
tor angles. In this same figure the corresponding ele- 
vator and tab angles are also plotted. In Fig. 8, 
curves of hinge moment coefficients for the elevator, 
with the tab floating in conjunction with the internal 
balance surface vented at the 10 per cent chord posi- 
tion, are shown along with the corresponding tab 
angles. In addition to these hinge moments, the pres- 
sure ratio P of the 10 per cent chord position is plotted 
against e for a, = 0. This pressure ratio agreed very 
well with the calculated pressures from Fig. 10 up to 
e = 15°, after which a drop in pressure gradient of about 
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Fic. 8. Comparison of measured and calculated elevator 
hinge moment coefficients and table angles for internal tab 
balance with 10 per cent pressure vents (R = —0.0426). 


50 per cent is to be expected. Fig. 8 also includes the 
calculated values of elevator hinge moments and tab 
angles as found from equations and charts presented in a 
later part of this paper. These calculations agreed 
very well except for the case of Cy, vs. e. In this case 
the calculations indicated a greater reduction in hinge 
moment than did the tests, when a, = 0. A possible 
explanation of this may be that the pressures created 
by the deflecting tab modified the pressures at the 10 
per cent chord position for the deflected elevator, the 
effect of which has not been included in the calcula- 
tions. 


METHOD OF PREDICTING HINGE MOMENTS FOR AN 
EXTERNAL ELEVATOR BALANCE 


The first problem in predicting the effect of an exter- 
nal elevator balance forming the forward portion of the 
stabilizer is to calculate the hinge moment coefficients 
for such a hinged surface. This has been done by 
graphically integrating the pressures of Fig. 10 to form 
Fig. 9, which represents the hinge moment coefficients 
of the forward portion of an airfoil having two serially 
hinged flaps. The hinge moment coefficients of the 
external elevator balance Cyy have been based upon the 
elevator chord rather than the balance chord; hence, 
these coefficients may be added directly to the hinge mo- 


ment coefficient for the elevator. The value of Cyy is 
given by the equation 


Cun = CuaCr + Cuno + Curee (29 
where 
Cr = ma + kp + ke) (30) 


For the present analysis the angle of attack a, of the 
horizontal tail is taken with respect to the fixed portion 
of the stabilizer and the angular movement of the ex- 
ternal elevator balance @y is measured positive for an 
upward movement of the balance nose. Here the ex- 
ternal elevator balance is assumed to transmit a frac- 
tion R of its hinge moment to the elevator which re- 
mains constant for all elevator angles; hence, the fol- 
lowing relations may be substituted in Eqs. (29) and 
(30). 


a= a, + by 
6= —by 
by = Re 


Cue = RCyw (due to stabilizer balance) 


If this external balance surface hinge moment is added 
to that produced by the elevator the resulting hinge 
moment of the elevator system becomes 


Cm = Crea: + Cet (31) 
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Fic. 9. Chart for calculating external elevator balance hinge 
moments. 
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where 
= Cuea + Cuam,R (32) 


a - Cree + [Cuskem, + Cupe|R -+- 
[Caml — ke) — Cype|R? (33) 


If my, is the slope of the tail lift coefficient curve for 
elevator free, then from Eqs. (30) and (31) for the condi- 
tion that Cy, = 0 


mi = mill + (1 — k)R+ RelCrea/Cuee (34) 


These equations permit the necessary horizontal tail 
characteristics to be calculated for an external elevator 
balance and are generally valid up toe = +15°. 

Calculated values of Cy,., and Cyee have been made 
for the model tested using the above equations and the 
results plotted in Fig. 5. These calculations are based 
upon Cy.g = —0.0031, Cyee = —0.0102, m, = 0.064, 
and k, = 0.59 as found from wind-tunnel tests. Then 
using Fig. 9 for E, = 0.30 and E, = 0.35, the value of 
Cus = — (—0.128) + (0.35)? = 1.045, Cupe = — 0.006, 
Cupe = —0.026, and k, = 0.91. 


METHOD OF PREDICTING HINGE MOMENTS FOR AN 
INTERNAL TAB BALANCE 


In this section equations will be developed permitting 
the elevator hinge moments to be calculated for the 
case where a hinged surface within the stabilizer is con- 
nected to float freely in conjunction with an elevator 
tab. 

The resultant air pressure acting on the internal sur- 
face for a one-degree change in a, or e may be found 
from Fig. 10 for any pressure duct location («/c). 
This figure has been calculated from reference (1) based 
upon characteristics for a sealed elevator. The pres- 
sure ratio difference acting on the internal surface 5S, is 
then 


P= Pia, + Pe (35) 


If c, is the distance of the center of pressure of this sur- 
face from its hinge and if R is the fraction of hinge mo- 
ment developed by this surface that is transmitted to 
the tab, the tab hinge moment coefficient due to the in- 
ternal surface will be 


Cus = (Paar + Pe) (Seo/S¢,)R (36) 


Here a positive value of R will develop a tab hinge mo- 
ment tending to deflect the tab to a positive value of 6. 
Combining the above equation with the normal tab 
hinge moment equation, the total tab hinge moment 
may be expressed as 


Cus = Cred: + Cusee + Cussd + 
(Paar + Pe)(Sp¢e/Sic,)R (37) 


The elevator hinge moment is 


Cue = Cres + Cree + Cresd (38) 
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Fic. 10. Chart for calculating chord-wise pressure difference 
ratio for an airfoil with deflected flap. 


By combining Eqs. 37 and 38 modified values of Cyeq 


and Cy,. may be calculated which incorporates the ef- 
fects of this floating tab. If these are denoted as 


Crea and Cy, then 
Cane ( SH) (39) 








a - Cava = CrreaC riba st 





CreaCuss CreaCuss\ Silt 
é.. = Cures E rie CuesCuse "= CuesPe (Ja (40) 
CueeCuss Cree "ome: | Sic 


The tab position in terms of a, and e may be written di- 
rectly from Eq. (37) for Cy; = 0 giving, 





Raat ea | Sas ‘ Z:(%#) R la = 
Cuss Cuss\ Sie 
Cue P, (*) | 
— + 00 \R | e (41) 
| = Cuss\ Sil: 


Knowing the tab position as a function of a, and e, the 
lift and hinge moments of the horizontal surface may 
be found for cases where the tab effectiveness reduces 
at the larger tab angles. This may be done by using 
wind-tunnel data obtained for various fixed tab set- 
tings. Calculated values of hinge moment coefficients 
and tab angles are shown in Fig. 8, based upon pres- 
sures from Fig. 10, and the following coefficients are 
determined from the wind-tunnel tests: 
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Cuea = —0.0048 Cys. = —0.0010 R = —0.0426 
Cuee = —0.0090 Cys. = —0.0010 k, = 0.59 
Cres = —0.0075 Cuss = —0.0097 mM, 0.073 


P, = 9.169 (from Fig. 10), P, = 0.058 (from Fig. 10) 


GENERAL DISCUSSION 


Several points concerning the application of this type 
of control will be discussed. At present there seems to 
be three desirable control force characteristics— 
namely, 

(a) Control forces should not exceed a given magni- 
tude. 

(b) Control forces should not reverse in direction, i.e., 
a pull should always be necessary to increase the air- 
plane angle of attack. 

(c) Motion of the control column should be positive, 
i.e., a pull on the control column should be accom- 
panied by a rearward movement of the control col- 
umn. 

Of these three, characteristics (a) and (b) seem im- 
portant; however, it is believed that (c) is not a vital 
necessity, since airplanes have been flown and reported 
satisfactory in cases where the control motion has been 
completely reversed, maintaining, however, positive 
control forces. In applying the present angle of attack 
balance, a reversal of control movement will occur 
when the elevator is forced to float in a positive direc- 
tion. One method of preventing this, if it is found 
necessary to do so, would be to provide a fore-and-aft 
horizontal tail, one surface being large enough for the 
pilot to handle with desirable control forces and the 
other forming a special stabilizer having an adjustable 
angle of attack balance. For this arrangement, as the 
airplane c.g. is shifted such as to reduce or increase 
the static stability, the stabilizing surface may be 
regulated to maintain the initial stability, giving the 
same control characteristics over an extended c.g. 
range. 

Another point of concern is the possibility of hunting 
of the airplane for elevator free when the airplane is un- 
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stable for elevator fixed.* This can only be studied by 
a careful dynamic analysis; but, if control friction 
forces are kept low and if a reasonable amount of angle 
of attack balance is always used along with proper 
damping from other parts of the airplane, this tendency 
should be avoided. 

A third point to be mentioned concerns the effect of 
compressibility shock waves effecting the airfoil pres- 
sures for high-speed airplanes. The external elevator 
balance may give trouble on high-speed airplanes be- 
cause of a change in pressure distribution and therefore 
will require a careful investigation. It should be re- 
membered that only 10 or 20 per cent of these pressure 
effects will be transmitted to the elevator; hence, con- 
trol forces may not be altered to any great extent. For 
the case of the internal tab balance it is believed that if 
sufficient damping action is provided, which may easily 
be done by restricting the area of the air vents, satis- 
factory results should be obtained at high speed. An- 
other combination, not tested but one that may be 
found desirable for a high-speed airplane, would be to 
connect the internal balance to the elevator directly in- 
stead of to the tab. This would require a greater in- 
ternal balance surface area, but for a reasonable amount 
of regulated control this would not be too large, pro- 
vided the elevator had some initial aerodynamic bal- 
ance. 

It should be pointed out that other control surfaces 
may likewise use this type of balance to an advantage in 
which case either a fixed or adjustable linkage may be 
found desirable. 
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* The following statement is quoted from page 33 of reference 
2: “Good results have been obtained with incidence stabilizing 
devices, which were tested on airplanes with very little or no sta- 
bility and which enabled the pilot to fly with release controls.” 
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Torsional and Aileron Flutter 


ZBIGNIEW KRZYWOBLOCKI* 


Polytechnic Institute of Brooklyn 


SUMMARY 

The aim of this paper is to find an approximate solution of the 
influence of various elements on the critical frequency in a binary 
system of torsional and aileron vibration. Assuming the system 
as a nonconservative one and taking into account the results ob- 
tained by von Karman and Biot, the influence of the kind of ma- 
terial and wing shape and airfoil thickness on the torsional fre- 
quency is considered. It is to be noted that damping is not 
taken into account. The idea of ‘“‘wing torsional rigidity” is in- 
troduced, being inversely proportional to the angle of twist at 
the end of the wing. The problem of influence of the kind of ma- 
terial, aspect ratio and span, wing taper and airfoil thickness on 
the torsional frequency is solved in a first approximation from the 
point of view of proportion. Finally, the influence of various 
materials of wing and aileron is taken into consideration. An 
attempt is made to solve this problem by the aid of the partial 
differential equation based on Newton’s Law. The equation is 
applied to an oscillating rod and permits the first natural fre- 
quency of the oscillating rod to be found. Using the Membrane 
Analogy, the differential equation is applied to a tubular mem- 
ber. The results are checked by the aid of Rayleigh’s Method. 
The following conclusions are drawn: 

(a) With regard to the influence of the kind of material on the 
binary system of torsional and aileron vibration, the decisive 
factor is the ratio of modulus of rigidity to the specific weight. 
High values of this ratio are obtained for steel and magnesium 
alloys; wood has the lowest ratio. 

(b) The influence of span and airfoil thickness are significant 
and the influence of wing taper is less important. 

(c) When a tubular member consists of various materials, a 
decisive factor is also the ratio of the modulus of rigidity to the 
specific weight. The following combinations might be advan- 
tageous: aluminum alloy or magnesium alloy wing and steel 
aileron, spruce-plywood wing and metal aileron. 

(d) Three different methods lead to the same result—that, 
from a system of torsional oscillation standpoint, the decisive 
factor is the ratio of the modulus of rigidity to the specific weight. 
The methods include the computations based on the magnitude of 
the torsional frequency given by von Karman and Biot, the dif- 
ferential equation based on Newton’s Law, and Rayleigh’s 
Method. 

(e) The results obtained cannot be applied directly to prac- 
tice because of the fact that damping, an important factor from 
wing flutter point of view, was not taken into account. 


SyMBOLS USED 


7 = kinetic energy 

K = Iw + ms? 

N = MgSWa 

R = Mete* 

le = moment of inertia of the wing about the 


elastic axis 
mass of the aileron 
= distance between the wing elastic axis and 
aileron hinge axis (the two axes are as- 
sumed to be parallel) 


Presented at the Aerodynamics Session, Eleventh Annual 


Meeting, I.Ae.S., New York, January 25-29, 1943. 
* Fellow, Polytechnic Institute of Brooklyn. 
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Wa 


la 
Maia* 


Pw 
Pa 
Fu P, 


Cu, Ci2, C21, Coz 


Rw, ka 


oe = 


Tmaz. 


e 


& YR 


O(x,t) 


Qs 


distance of the center of gravity of aileron 
from the hinge axis (positive if the center 
of gravity is behind the hinge axis) 

radius of gyration of aileron 

moment of inertia of aileron about the 
hinge axis 

angular deflection of wing relative to di- 
rection of flight 

angular deflection of aileron relative to di- 
rection of flight 

generalized forces 

coefficients 

coefficients 

torsional frequency of oscillation 

torsional frequency of the system, when 
aileron mass is concentrated on the hinge 
axis 

frequency of aileron oscillating around 
hinge 

flight speed 

angle of twist at end of wing 

torque 

area enclosed by the mean curve of the 
periphery 

modulus of elasticity in shear (modulus of 
rigidity) 


= Poisson’s ratio 


i] 


wall thickness 

length of centerline of ring section 

length of wing (semispan) 

shearing stress 

ultimate shearing stress 

wing torsional rigidity 

specific weight 

aspect ratio = L/; = F/, = L/F? 

bottom side of ‘‘wing working section’”’ hav- 
ing rectangular shape 

height of ‘‘wing working section”’ 

length of wing chord 

wing area (semispan) 

radius of gyration of wing 

weight of wing 

weight of airplane 

constants 

coefficient depending on kind of flight 

constant 

design load 

length of wing root chord 

length of wing tip chord 

t,/te = taper ratio 

length of intermediate chord 

a/b—assumed to be inversely proportional 
to airfoil thickness 

thickness of wing airfoil 


= angle of twist of section at distance x from 


the origin relative to the section at the 
origin 

twist per unit length 

constant depending on the shape and size 
of section of rod 
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I = moment of inertia of cross section of rod 
about the x-axis (centroidal moment of 
inertia) 

g = acceleration due to gravity 

t = time 

p = (y/g)I = moment of inertia per inch length 

C1, C2, Cz, C4 = constants 

w) = the first natural frequency 

l = length of rod 

U = potential energy 

Oo = angle of twist at free end of rod 


INTRODUCTION 


_ A POSSIBLE COMBINATION of flexural-torsional 
and aileron oscillation probably the most critical 
is the case of torsional and aileron oscillation. Strictly 
speaking, the influence of wing and wing shape on os- 
cillation refers more to the problem of pure torsional or 
flexural oscillation than to the case of binary system of 
torsional and aileron oscillation. 

In the latter case an important factor is the influence 
of aerodynamic qualities of the airfoil and the elasticity 
of the control system. In practice there probably 
never is a pure case of torsional or flexural oscillation 
but chiefly a general case of flexural-torsional and 
aileron oscillation. For this reason the author gives 
main consideration to the case of the binary system of 
torsional and aileron oscillation. 

THE INFLUENCE OF CONSTRUCTIONAL MATERIAL AND 
SHAPE OF WING 


Basic Assumptions 


In the following the wing is represented by a rigid 
cylindrical body, having a streamline section that is 
elastically restrained to the wing elastic axis. The wing 
is clamped to the airplane fuselage, which is considered 
a rigid base. The aileron is hinged at the aileron axis. 
Special attention is paid to the case of torsional oscilla- 
tion of wing combined with oscillation of the aileron, 
because, as is known, this vibration is usually more ef- 
fective in its results on account of the change of angle 
of attack. 

The kinetic energy of the wing-aileron system can 
be expressed in the form: 


T = (Kpu? + 2Npuba + Rba?)/2 (1) 


Assuming the system as a nonconservative one, the 
moments as the generalized forces, and the angular 
deflections as the generalized coordinates 


P, = —CupPw — Ci2pa (2) 


P, = —Crupy — Copa (3) 


Cup. represents the restoring moment about the 
wing elastic axis produced by an angular deflection of 
the wing. Its magnitude depends on the torsional stiff- 
ness of the wing and on the aerodynamic moment 
caused by a change of angle of attack. 


AUTICAL 
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Cp, represents the moment about the wing elastic 
axis caused by the aileron deflection. 

CaP.» represents the hinge moment caused by the 
change of wing's angle of attack. 

CP, represents the hinge moment caused by the 
change of angle of attack of the aileron. 

The system is nonconservative because: 


Cy F C2 (4) 


Substituting the above equations in Lagrange’s equa- 
tions 


Kbw + Nps = — CuPw — Ci2Pa (5) 

Nw + Rea = —CaPo — CoxPa (6) 
Substituting 

Pu = hyd (7) 

Pa = kei (8) 

wy? = Cy/K (9) 

Woo” = C22/R (10) 

and assuming that w, = O—i.e., that the center of 


gravity of the aileron lies on the hinge axis, that the 
aileron is statically balanced, and that there is no dy- 
namic coupling between the coordinates p, and p,— 
von Karman and Biot! obtained the following condi- 
tion for the stability of the system: 


2=— on” + waa" si 


2 


~ 


y(2" w ) — w117029” + ts (11) 


The system is stable only if w is real—i.e., w* is real and 


positive. To accomplish this there are two conditions: 
on? = We29" 
fou? + wm)? gg | Cnn (49 
(2 2 oi Wi woo” of KR ( ) 
and 
. A . =, Sam ’ 
(sot on’ at = ) — w1"weo2" + a > 0 (13) 
or 
(wn? ina We9”)? + 4 (C122 /KR) > 0 (14) 


Usually the product of cy. and cq is less than zero. 
w11” should be as large as possible. The system is un- 
stable in the neighborhood of w;; = wy. The case wz + 
0 is not taken into consideration because the results ob- 
tained will not be changed in that case. Both w, and 
w2 depend on the flight speed; the coefficient c,, and 
the torsional frequency w, decrease with the flight 
speed; the frequency w2. increases with the speed. 
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At a certain speed these two frequencies will, of course, 
coincide. 


Wing-Torsional Rigidity 


All four coefficients cy, C12, C21, C22 given above depend 
on the aerodynamic properties of the wing and the 
aileron. Besides this, only one coefficient, ¢,, depends 
on the constructional properties of the wing. To con- 
sider the influence of the wing material and wing shape 
on the torsional frequency w,, the influence of the above- 
mentioned elements on this frequency must be taken 
into account. To accomplish this and to avoid compli- 
cated computations, the problem will be considered from 
a qualitative and not from a quantitative point of 
view. To do this, the conception of proportionality 
instead of the conception of equality is intro- 
duced. 

Torsional rigidity is usually the factor by which the 
torque is divided to obtain the twist per unit length. 
In these considerations another element introduced 
concerns the torsional rigidity of an airplane wing. 
In wing construction the most important consideration 
is perhaps the angle of twist at the end of the wing. 
Wing twisting causes a change of aerodynamic quali- 
ties. In many countries the aircraft requirements often 
limit this angle to a certain value, mostly 2°. As the 
measure of the ‘‘wing torsional rigidity” unity divided 
by the greatness of this angle is accepted. Using 
Bredt’s equation for a thin tubular section 





o- [ ie © al (15) 
This equation can be modified to 
M, = 2Aér (16) 
or 
M, = 2Abtmaz, (17) 


Substituting Eq. (17) in Eq. (15) and assuming, for 
comparison, 6 = a constant: 


L 
‘ 
6 = — ds dl 18 
Lat as 
As the wing torsional rigidity, it is assumed that 


= 1/0 (19) 





The Influence of the Kind of Material 


It is assumed that: (a) the aerodynamic properties 
of the wing are the same (A, airfoil, and chord); (b) 
the shape and construction of wing are the same; 
(c) the change of material does not affect the values 
A, s, and the radius of gyration of the wing; and (d) 
damping is not taken into account. 

Using these assumptions 





(20) 


TORSIONAL AND AILERON FLUTTER 





From Eq. (9) 


P Cu Cn 
Se 21 
. K  Iy + m,s* @n) 





It is known that the restoring moment cp, corre- 
sponds to the torsional rigidity of the wing and to the 
aerodynamic moment caused by a change of angle of 
attack. In the first approximation, not considering the 
degree of proportionality, only the torsional rigidity 
and weight of the wing can be taken into account. 

@,? = const. .. = euteitices ~ =. (22) 

Le W. WwW. 

The weight of the wing is proportional to 6 and y. 
From Eq. (17) the assumption that M, = const. and 
A = const. gives 


5 = const. = 2 ~ = (23) 
T maz. T mar. 
S@bstituting Eqs. (20) and (23) in Eq. (22) gives: 
@117 = const. G ~ G (24) 
Y Y 


Therefore, in considering the influence of the kind of 
material on the binary system of torsional and aileron 
oscillation, the decisive factor is the ratio of modulus of 
rigidity to the specific weight. The values of this ratio, 
given in Table 1, are taken from references 2 and 6. 





Taste 1 
Ultimate 

Shear 

Modulus of Specific Stress, 

Rigidity,G, Weight, © ) Tmas- 
Kind of Lbs. per y, Lbs. Y Lbs. per G/- 
Material Sq.In. per Cu.In. 10,000 Sq.In. Tmas- 
Spruce 84,000 0.0157 536 1,000 84 

Spruce-ply wood 

45° 420,000 0.0188 2,235 2,370 177 


Mild 1025 steel 11,000,000 0.2833 3,890 35,000 314 


C.M. steel 11,000,000 0.2833 3,890 55,000 200 
17ST aluminum- 3,800,000 0.1010 3,760 33,000 115 
Alloy sheet 
Magnesium al- 
loys* 2,500,000 0.0660 3,790 18,000 139 
Deaeeteae 8 8}}§©=«—COswu ees ionn 3,380 sien he 





*According to reference 6 the modulus of elasticity of magnesium alloys 
is 6,500,000 Ibs. per sq.in. Assuming G = E/2(1 +1) andr = 0.3 the above 
given value is obtained. To obtain rmaz. for magnesium alloys it is assumed 
that 

Tmax. = 0.6 X ultimate tensile stress 
This ratio in a first approximation is in agreement for the metals used; ac- 
cording to references 6 and 7 the ultimate tensile stress for magnesium 
alloy structural shapes equals 30,000 Ibs. per sq.in. 

¢ According to reference 12 tests for Duramold show: 

3,060 = E/1,000 X specific gravity 


and as 
Specific gravity = 1,728 X specific weight/62.5 
Assuming 
G = E/2(1 + r),r = 0.25 
one has 


G 3,060,000 X 1,728 
62.5 X 2.5 





= 3,380 X 10 


specific weight 





The value obtained presents only the first approximation and the results 
should perhaps be used cautiously. 
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As is seen, metals have a greater ratio of G/y in 
comparison with wood, even plywood. The high value 
of magnesium alloys is interesting, as is that of Dura- 


mold. 
If, in a special case, it is assumed that the weight of 


the wing is constant Eqs. (20) and (22) give 
G 


T maz. 





wi” = const. 8 ~ (25) 
Eq. (25) is in agreement with Eq. (4: 70) given by 
Younger.? Some values of this ratio are given in Table 
1. Again the great advantage of metal is seen, especi- 
ally steel, in comparison to ordinary wood. The high 
value obtained for plywood is striking. 


The Influence of Aspect Ratio and Span 

Assume that (a) the airfoil section and wing surface 
area are constants; (b) the shape of the wing is rec- 
tangular; (c) the geometric details of construction are 
the same regardless of the materials; (d) the spar or 
“working section” of the wing is a rectangular shell 
having lower side a, height }, and a uniform thickness 
of walls 5; (e) the side a is proportional to length of 
chord; (f) the height } is proportional to thickness of 
profile and also to length of chord (thickness is ac- 
cepted as a percentage of chord length). 

For comparable aims Eq. (18) may be modified to 
take into account an intermediate section, approxi- 
mately in the center of the wing span. This assumption 
and the equalities 


fds =s22(a+ 5) (26) 
A =ab (27) 


I|S 
S 


permit the writing of 


e = T maz. ue b)L (28) 


where a and 3 refer to the accepted intermediate sec- 
tion. As 


a = const.t ~¢ (29) 
b = const.t ~t 
therefore 

(a+ 6) = const.t ~t (30) 
ab = const. t? ~ ?? (31) 

2 
ee a ene ee (32) 

t t F 


& =const. 1/A = const./L?~1/L* — (33) 
Modifying Eq. (21) 
a 7 ea (22) 
Ty + ms? Ie Ie 


the change of J,, must be taken into account. Assum- 
ing that the distance measured along the wing chord 


, 
ay = 
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line is proportional to the length of the chord line gives 
Zn? = const. t? ~ (1/L*) (34) 
There are various formulas expressing the change in 


wing weight with change of span. For instance, in 
reference 3 





‘a i W ue 
- = (1, + ko ) nL F 4 kz (35) 
In reference 4 
W, 92 =. (36) 
1,000¢ 
or 
W,=K ee (37) 
1,000¢* 


It is seen that as a first approximation® 


W,, = const. L? ~ L? ~ - (38) 
Therefore 
Ip ~ We X ig’ ~L~= (39) 


Substituting Eqs. (33) and (39) in Eq. (22): 


wy" ~~ ® ~ i ~~ ih (40) 
i 
It can be seen that the influence of span on the binary 
system of torsional and aileron vibration is significant. 
The author obtained similar results by a graphic 
method.°® 


The Influence of Wing Taper 


The assumptions are: (a) the aerodynamic profile 
and area of wing are the same; (b) the construction and 
material of wing are the same; (c) the plan form of the 
wing is a trapezoid; (d) the length of the wing root 
chord ¢; is constant; other assumptions are similar to 
those applied above. 

It follows that 


ty = ne (41) 
F = Lt,, = constant (42) 
— ti + te a (1 + (1/n)] (43) 
2 2 
1 1 


(44) 





L = F/t, = const. ~ 
1+ (1/n) 1+ (1/n) 


Substituting Eq. (44) in Eq. (40) gives 
wy? = const. [1 + (1/n)]* ~ [1 + (1/n)]® (45) 
The function 
L = fin) (46) 











(36) 


(37) 


(38) 


(39) 


(40) 


nary 
‘ant. 
phic 


‘ofile 

and 
f the 
root 
ir to 


(45) 


(46) 
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has a horizontal asymptote. It is seen that the in- 
fluence of taper is greater for smaller values of m and 
less when m increases. When two wings are com- 
pared, one rectangular and the other trapezoidal, hav- 
ing the same wing area, the results show that (a) 
assuming the same value of A and changing /; and fs, 
both wings are equivalent from the standpoint of 1; 
(b) assuming the same length of wing root chord, the 
trapezoidal wing is less advantageous than the rec- 
tangular one (because J is greater). 


The Influence of Airfoil Thickness 


The assumption is made that the only dimensional 
variable is the thickness of airfoil and that the height 
b of the working section is proportional to the thickness. 
Modifying Eq. (28): 


0 = const. ng const. oll + (1/x)] ~util 
ab a*(1/n) 
(47) 


Therefore, since @ is always large in comparison 


to b 
Om u (48) 


Assuming that the weight of wing is affected by the 
change of airfoil thickness in the manner given by 
Younger’? [Eq. (4: 40)]—.e., that the wing weight is 
inversely proportional to the thickness of the wing, 

1 1 

W. = const. — ~ - (49) 

H H 
uw is inversely proportional to the airfoil thickness. 
Substituting Eqs. (19), (48), and (49) in Eq. (22) 


T 
wy? ~w a ~ H? (50) 
M 
A significant influence of airfoil thickness is seen es- 
pecially when taking into consideration the tendency 
to use thin airfoils on modern airplanes. 


THE INFLUENCE OF VARIOUS MATERIALS ON WING 
AND AILERON 


In practice, the wing frame sometimes consists of 
various miaterials—for instance, dural and _ steel. 
Also, different materials are used for wing and aileron, 
stabilizer and elevator, fin and rudder. 

In a binary system of torsional and aileron oscilla- 
tion, probably the deflections of the ailerons around 
the hinge line as a whole are more important than the 
torsional deflection of the aileron in itself. But one 
can imagine that the aileron is not free and that the 
elastic restraints such as control wires are not flexible 
enough. On the assumption that one applies the ir- 
reversibility of aileron controls, the wing-aileron system 
can be treated as one unit. In the design office, the 
Opinion is sometimes held that the same material for 
both parts increases the possibility of wing flutter. 
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In other cases the possibility of wing flutter is more 
remote. In the case of the same material the parts are 
“in phase”’ and in another case “‘out of phase’ from the 
wing oscillation standpoint. Evidently the “‘critical 
torsional frequency”’ of the system is removed in the 
second case to a higher value. 

In an attempt to solve this problem the torsional 
oscillation of a rod fixed at one end and free at the 
other will first be considered. Newton’s Law will be 
applied to obtain partial differential equation of the 
motion. Damping is not taken into account. 

Assuming that the x-axis is taken along the centers 
of inertia of the cross sections of the rod, that there is 
pure torsion in each element of the rod, and that the 
section at distance x from the origin is twisted through 
the angle 0(x,/) relative to the section at the origin, 
where w denotes the angular velocity of this section, 
an element of length 6x is twisted through the angle 
60(x,t). The twist per unit length of this element is 


y = 00(x,t)/Ox (51) 
If M, denotes the torque in this element 
M, = CGy (52) 


The constant C depends on the shape and size of the 
section of rod and possesses the dimensions of the 
moment of inertia of an area. Assuming that J de- 
notes the moment of inertia of the cross section of rod 
about the x-axis, the moment of inertia of the element 
of the rod of length 6x is 


(y/g)L6x (53) 
The angular momentum of this element is 
(y/g)Iwix (54) 


The rate of increase of the angular momentum is equal 
to the total torque on the element. The total torque 
is the excess of the torque at (x + 6x) over that at x, 
also: 





(OM,/dx) dx (55) 
and 
© ((y/eylwtx) = Me 5, (56) 
ot . 6x 
dw oM 
I — 6x = ; 57 
(y/g) oa seer ae (57) 


Taking into account Eqs. (51) and (52 


O°O(% 1) oy DOCH) 


58 
ot? Ox? (68) 


This equation is in accordance with methods given in 
references 8 and 9. The variable in this case is the 
angle of twist O(x,/). Assuming that the rod vibrates 
harmonically at some natural frequency 


Q(x,t) = O(x) sin wt (59) 
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. i zi (52 
o°(@.1) ee ' < eeltaia at! (60) Therefore the expression for CG in Eq. (52) becomes 
or" } S (4A°G/ds)é (75) 
O*O(x,t) _ d*O(x) sin wt (61) However, the thickness 6 is inversely proportional to 
Ox? dx? the maximum allowable stress Tm, Thus, Eq. (74) 
Substituting Eqs. (60) and (61) in (58) becomes 
. Z dO(x : ae i A le 76 
pO(x)w? + CG = = 0 (62) M, aw f * ie tee: (76) 
« Assuming that the dimensions of the wing are un- 
changed 
ee? 4. ais) = 0 (63) G 
dx? CG M, = cw - (77) 
T maz. 


This is an ordinary differential equation of the second 
The general solution is 


O(x) = co sinx \"" + ccosx y=" 64 
se" _ 


degree. 


Assuming that the function 





O(x) = f(x) (65) 
is a quarter sine wave leads to the end conditions 
x = 0, A(x) = 0 (66) 
x = 1, dO(x)/dx = 0 (67) 
From the first condition, 
cf = 0 (68) 
O(x) = aq sinx pow? (69) 
CG 
The second condition steed 
dO(x 2 2 
-_ =o We cos! pa* = 0 (70) 
dx CG CG 
or 
nye! = 2/2, 38/2, 5e/2.... 71 
CG /2, 3x/2, 5a/ (71) 


The first natural frequency 


wo, = /2 y¥& (72) 
l p 


Applying this to a thin tube or to the case when the 
cross section of a tubular member has more than two 
boundaries, by use of Eqs. (146) and (162) in reference 
11, one can write, on the basis of the Membrane Anal- 
ogy: 

4A°G 


5 (73 
— (73) 


M, = 





It is seen that in a first approximation Eq. (52) can 


be modified to 
M,=y 7= a °C (74) 








Let it be assumed that the tubular member consists 
of various materials, ideally connected together. In 
computing the first natural frequency of this system, 
Eq. (72) can be applied substituting instead of CG the 


expression 
G Gy aie 
af = = & eae (78) 
T mar. = T \maz.i 


i= 





where the subscript 7 refers to various materials. 
Similarly, Eq. (53) can be modified. So instead of 


(y/g) I 6x (53) 
the following expression may be introduced: 


[(na/e)li + (2/g)la +... .(%9/¢Malox = 
| = (eh fi (79) 


if y1, Yo, .--- Yn, refer to the first material, second ma- 
terial, and so on, and J, J, .... I, signify the partial 
moment of inertia of the parts of the tubular member 
made from the first and second material, etc. The 
value J, in a tubular member is proportional to the 
corresponding wall thickness 6, which is in turn in- 
versely proportional to Tar, Therefore 


i=n 


) (y:/g)I; = const. Ye (80) 
T maz.; 
i=l oe 


+= 


In computing the first natural frequency of torsional 
oscillation of a tubular member consisting of various 
materials, one must use instead of Eq. (72) the equa- 
tion: 


oy, = — an’ = 
Pee i : » Pies i 


Y nieth: (81) 


In practice a wing and an aileron do not represent a 
rigidily connected system as the one assumed in the 
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above given computation, because the wing-aileron 
system is an elastically connected one. However, the 
above given computations give an approximate idea 
with regard to the “‘critical torsional frequency” of a 
system consisting of various materials. 

It can be seen that the frequency « is influenced by 
a selection of materials having different ratios of G/y. 
Table 1 shows that such combinations as plywood and 
steel or aluminum alloys and steel might be advan- 
tageous. In practice, these combinations are applied. 
The result, obtained above, can be verified by Ray- 
leigh’s Method. In any free vibration the following 
condition must be realized :"° 


(Total kinetic energy) + (Total potential energy) = 
Const. (independent of #) (82) 


In other words, this condition may be expressed in the 
following: 


maximum kinetic energy 
(83) 


maximum potential energy = 


The energies of the oscillation rod (torsional oscilla- 
tion), having circular cross section will be computed. 

The total potential energy stored in a slide dx of the 
rod can be expressed: 


dU = M,dx/2GI (84) 
From Eq. (52) 
M, = 1G 24) (52a) 
ox 


The maximum potential energy occurs in the most de- 
flected position. Taking that moment into considera- 
tion, it is seen that 0 will be independent of ¢ (time). 
One has to assume a shape for the first normal deflec- 
tion curve, i.e., angle 0 versus distance x. As before, a 
quarter sine or cosine wave may be assumed. Choosing 











O(x,t) = O(x) = Op sin a (85) 
remembering that 
OO!) _ PC) _ oF cos (86) 
Ox dx 21 21 
and substituting Eqs. (52a) and (83) in Eq. (87) 
lw 2 
U = [ es A (87) 
0° 26r 
gives 
U = “5 GI®,? i (88) 


The maximum kinetic energy occurs in the undeformed 
position where the velocity is greatest. The kinetic 
energy of a rod element dx is: 


dT = (p/2) 6°(x,t) dx (89) 
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But 
8 = w0(x,!) (90) 


Substituting Eq. (85) into Eq. (89) and integrating 


Tt = {1% pe? 042(sin = dx (91) 
0 2 21 
After integration 
T = plw®,2/4 (92) 
From the condition 
T=U (93) 
one has 
w ee (94) 


Eq. (94) is identical with Eq. (72). 

The assumed function does not fulfill satisfactorily 
the required end conditions. The assumed curve 
must be horizontal at x = 0 and the angle of twist at 
the free end must be equal 6). The following equation 
satisfies these conditions: 


O(x) = Op (1 — cos =) 


= (85a) 


This function does not present in the present case an 
approximation of a high degree. The result is: 


1 
U = — GI0,?x? 88 
 paeeaeaaia ne 
= [(?/4) — (?/1)]plw*@,? (91a) 
0.74 1G 
1 ——— ye (94a) 
l p 


Another curve may be assumed, but since this problem 
has been treated from a general point of view, this will 
not be discussed any further. 

When applying the above method to a thin tube or 
to a tubular member having more than two boundaries, 
Eq. (73) instead of (52a) should be substituted. The 
same result would be obtained as formerly. Of course, 
in any given case it is necessary to take into account 
the change of J and p along the span. This will have 
an influence on the value of U and T. 


CONCLUSIONS 


The following conclusions may be noted: 

(1) As far as the influence of the material on the 
binary system of torsional and aileron vibration, neglec- 
ting damping, is concerned, the decisive factor is the 
ratio of modulus of rigidity to the specific weight. 

(2) In connection with (1) the following materials 
might be used, listed in the order from the most to the 
least advantageous: steel, magnesium alloys, 17ST 
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aluminum-alloy sheet, Duramold, spruce-plywood 45°, 
spruce. 

(3) In a special case, assuming that the weight of 
wing of various materials is constant, the decisive 
factor in the binary system of torsional and aileron vi- 
bration is the ratio of modulus of rigidity to the ulti- 
mate shearing stress. 

(4) In connection with (3) the materials might be 
listed in the following order: 1025 steel, C.M. steel, 
spruce-plywood 45°, magnesium alloys, 17ST alumi- 
num-alloy sheet, spruce. 

(5) Worthy of note are the high values obtained for 
steel in comparison to other materials. 

(6) From (2) it can be seen that magnesium alloys 
might be advantageous and that plywood is less ad- 
vantageous. Duramold is better than plywood. 

(7) 17ST aluminum-alloy sheet has an average value 
from a torsional aileron vibration standpoint. 

(8) The influence of span is significant, the square 
of the torsional frequency of the system being in- 
versely proportional to the cube of the span. 

(9) The influence of taper is greater for smaller 
values of taper ratio and less for greater values. 

(10) The influence of airfoil thickness is also signifi- 
cant, the torsional frequency being proportional to the 
thickness. 

(11) When a tubular member consists of various 
materials, the decisive factor is again the ratio of modu- 
lus of rigidity to the specific weight. (Damping is not 
taken into account.) 

(12) In connection with (11) the materials might be 
listed as follows: steel, magnesium alloy, 17ST alumi- 
num-alloy sheet, Duramold, spruce-plywood 45°, spruce. 

(13) From the above it may be seen that the follow- 
ing combination might be advantageous: 
ered by 17ST aluminum sheet or by magnesium-alloy 
sheet—-C.M. steel or 1025 steel aileron; wing covered 
by spruce-plywood 45°—C.M. steel, 1025 steel or 
17ST aluminum alloy sheet aileron. 
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(14) From standpoint of torsional vibration the de- 
cisive factor is the ratio of the modulus of rigidity to the 
specific weight when damping is not taken into ac- 
count. 

(15) It is characteristic that the last conclusion is 
based on three different methods: i.e., the computa- 
tion of the torsional frequency given by von Karman 
and Biot, by Newton’s Law, and by Rayleigh’s Method. 
All lead to the same result. 

(16) The above results cannot be applied directly 
in practice, because damping, an important factor from 
wing flutter point of view, was not taken into account. 
Damping probably would change the above results 
considerably. 
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